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Chapter 1 

Introduction 



In 2008, the Large Hadron Collider (LHC) at the Conseil Europeen pour la 
Recherche Nucleaire (CERN) will start operating. The LHC is the next step in 
discovering the mysteries of particle physics and the fundamental forces. Reaching 
energies at the level of 14 TeV, the discovery of physics beyond the SM and new 
particles will become possible. Since the Standard Model (SM) of particle physics 
is most probably not the final theory, searches for theoretically motivated models 
have already begun and physicists wait for their confirmation. 

But making sure that one really observes a specific model is far from trivial. 
Two very promising candidates for physics beyond the Standard Model are Super- 
symmetry (SUSY) ill] and Universal Extra Dimensions (UED) [2,]- An important 
difference between them is their difference in the spin of the particle fields. In 
Supersymmetry, the usual space-time symmetry of the SM is enlarged towards a 
more general supersymmetry algebra by adding operators that change the spin 
of the particles. There are theoretical bounds like the Higgs mass, which allow 
to exclude or confirm SUSY at the LHC to a high certainty level. The model of 
Universal Extra Dimensions is a string theory inspired, extended version of the 
well known Standard model. Allowing all particles to propagate in five or more 
dimensions produces new particles. In the theory they appear after performing 
compactification of the extra dimensions. 

To show that SUSY is indeed the theory which is realized by nature, one has 
to measure the spin of all new particles. Since SUSY particles are not stable, 
except for the lightest supersymmetric particle, they are treated as parts of long 
decay chains. In this thesis we concentrate on the spin of the gluino, which is a 
Major ana fermion in SUSY. Major ana fermions are interesting anyway, since they 
can produce like-sign dilepton signatures which are a clear sign for new physics, 
as explained in Using like-sign dileptons in searches for new physics yields 
a much clearer signature than the one from missing energy measurements. But 
since like-sign dileptons can also be produced from the decay of a boson with an 
adjoint color charge, measuring the spin is crucial. What has to be done is to 
compare SUSY to a model with a "bosonic gluino", assuming equal masses of the 
particles. We therefore consider UED, where the gluino is substituted by a heavy 
Kaluza-Klein-gluon (KK-gluon) with bosonic spin-statistics. 

In recent publications the topic of spin measurements in decay chains gained 
higher popularity. In [3], it is shown that measuring the spin of the squark is 
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possible by comparing angular and charge asymmetries of the decay products in 
UED and SUSY. 

It is demonstrated in [B], that the spin of the gluino has a measurable impact 
on kinematic distributions and angular correlations of the decay products in gluon- 
gluon and gluon-quark collisions. Therefore leptonic and hadronic correlations of 
the decay products are investigated in UED and SUSY. It is presumed that the 
difference in the boost of the gluino and the KK-gluon is the main reason for the 
differences in kinematics of the decay products. As a second influence the coupling 
structure is mentioned. 

In 0], it is studied if UED can perfectly imitate the gluino decay, not necessarily 
assuming equal masses for the particle spectrum of both theories. It is shown for 
the decay of the gluino into two jets and a neutralino, that UED and SUSY can be 
differentiated, assuming that the gluino is lighter than all squarks. The influence 
of the coupling structure is investigated in a toy model and a comparison of the 
invariant masses for SUSY and UED is performed. The number of events needed 
to discriminate both theories is calculated. 

In this thesis we investigate the process of gluon-quark collision, producing 
two 6-jets, a light quark jet and two neutralinos in SUSY. An equivalent process 
is calculated for UED and all cross sections are normalized to make a comparison 
possible. Kinematic distributions of the decay products are compared in SUSY 
and UED, aiming to explain the difference in the angular distributions between 
the outgoing 6-jets. We study whether the influence from the different boost of 
the gluino and KK-gluon or the different coupling structure is more relevant to 
the correlations of the final state bottom jets. Therefore we vary the Kaluza-Klein 
mixing angle and numerically eliminate the effect from the difference in the boost 
distributions. In comparison to [6], we assume that the gluino is heavier than all 
squarks as it is e.g. in the SPS la scenario. 

This thesis is structured as follows: In the second chapter, we review the 
definition of the cross section and sum up a few important concepts. 

In the third chapter, we shortly introduce SUSY, by giving an impression of 
how a supersymmetric Lagrangian in general looks like in the spinor notation. 
We do not extensively review the complete particle spectrum of SUSY here. The 
theory of UED is introduced in somewhat more detail and all Feynman rules and 
Lagrangians are calculated completely, since literature on this is more rare. 

In the fourth chapter, we present our results for squark and gluino production 
cross-sections. The partonic cross sections are calculated analytically for all SUSY- 
QCD processes. Hadronic cross sections and kinematic distributions for the LHC 
are presented numerically for the case of gluino-squark pair production. The 
partonic cross-section for the production of a Kaluza-Klein gluon-quark pair in 
UED is also calculated analytically. This process is compared to the corresponding 
SUSY process at the production threshold. 

In the fifth chapter, the general kinematics of decay chains are discussed and 
the concept of our program is explained. We explain how the decay phase space 
is tested. 

In the sixth chapter, calculations are performed for the case of a SUSY and 
UED decay chain. We investigate the 2 — > 5 process of quark-gluon collision in the 
two models and compare them. Since mass spectra are equal and cross sections are 
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normalized, their only differences come from angular correlations, deriving from 
the spin of the particles. Results are compared for SUSY and UED. The influence 
of a varying mixing angle in the UED couplings is investigated. We also discuss 
the origin of the differences in the angular distributions which were claimed to be 
due to the different boosts of KK-gluon and gluino or to the different coupling 
structure of the heavy quark partners. 



5 



Chapter 2 

Basic Concepts 



2.1 Calculation of Cross Sections 

The total cross-section is an important theoretical result of high-energy particle 
physics since it gives an observable quantity which can be compared to experi- 
ments. Differential distributions of the cross section can equally well be calculated 
and observed at colliders. 



2.1.1 Definition of the Cross Section 

In general the differential cross section is defined as the product of transition rate 
per scattering center and the number of final states reachable, divided by the 
incident flux. The transition rate from an initial state \i) to a final state (/| is 
given by 

where T is a finite time interval and Sfi is called the S-matrix element. In a 
process where two initial particles with momenta pf = {Ei,pi), i = 1,2 collide 
and N particles with the momenta p'^ = {E'^^p'^j, f = 1,. . . ,N come out, the 
matrix element Ai is connected to the S-matrix element by 



To get from these relations to the transition rate to final states with momenta in 
the intervals {p'^p'^ + dp'^), one has to multiply by 

T-r Vd^p'f 
/ ^ ^ 

where the volume V is considered to contain exactly one scattering center. The 
incoming flux is Vrei/V, where Vrei is the relative velocity of colliding particles. 
Combining everything one finds 
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2.1.2 Parametrization for 2 — > 2 Processes 

For kinematical analysis it will be important to find the differential distributions 
of cross sections in variables convenient for observation. Therefore distributions 
of transverse momentum pt and rapidity y will be used. We show how they can 
be derived in the case of two massless colliding particles with momenta ki,k2 
producing two massive particles with momenta pi , p2 as in the case of SUSY. The 
usual Mandelstam variables are given by 

S=ikl+k2f, t={k2-p2?, U={ki-P2f. (2.5) 

In 0], the differential cross section is given as 

da 1 



dt 167rs2 



\M{s,t)\'. (2.6) 



When the beam axis is parallel to the z-axis, the variables pt = {p^. + Py)^^'^ and 
y for an outgoing particle with momentum p2 and mass m2 are given by 



P? = m? - ml] sin^r , (2.7) 

(2.8) 



2 



P2 + P2,L 
P2 - P2,L 



where 6* and the longitudinal momentum pi = \p\ cos9* are defined in the center 
of momentum system (CMS). A very comfortable feature of pt for practical cal- 
culations is its independence of boosts parallel to the beam axis. The advantage 
of the rapidity is its additivity under boosts. In the CMS energy and momentum 
of the outgoing particles are given by 

s + mf-ml ^ ^ X^{s,ml,ml) 

= , p* =p* = — ^- — 11 — ^ (2.9) 

^ 2^5 ' yi ^^2 2^5 

where 

X{x, y, z) = {x-{^+ V^f) (x - ( - V^f) (2.10) 

and E2 is obtained by exchanging the masses of the outgoing particles. For mass- 
less incoming and massive outgoing particles with masses mi and 1712 one finds 

t = ml-^{s + mi-ml) + ^{s + ml-m^^)(^^^\ , (2.11) 



which leads to 

dt (s + — mf) 



dy 2 cosh^T/ 



(2.12) 



Finally the differential rapidity distribution for the particle with momentum p2 is 
given by 
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To derive the differential pt distribution we use 



t = m\ — -(s + m\ — m\) + - A2 (s, m^, m|) cos6* , 
which is easily derived in the CMS, and with ()2.7p we find 



ti/2 = m\ - ^(s + - ms) ± ^Aa (s, ml,ml) 



1 



(2.14) 



(2.15) 



where ti corresponds to the case of the "+" sign with < < ^ and ^2 corre- 
sponds to the "-" sign with ^ < < vr. The resulting differential cross section of 
the transverse momentum is then given by 



dcT 
dpt 



dt 



dpt 



da 



1 



Pt 



dt 327rs2 ^A(s,m^,m^) 



Pt 



\M{s,t = h)\'' + \M{s,t = t2)\' 



Later we use histograms to calculate differential cross sections for any variable 
within just one run of our program. A parametrization for differential cross sec- 
tions of hadronic 2^2 processes is given in section 14.31 An overview of particle 
kinematics can be found in [3]. 



2.2 The Parton Model 



In this section the parton model is introduced. The parton model links the cross- 
sections on the hadronic level to the collisions of the constituents of the hadrons, 
quarks and gluons, also called partons. 

This is a fundamental task for making predictions for experiments at hadron 
colliders. Effectively, the parton model describes the collision of hadrons by the 
collisions of single partons while the so-called spectator partons are not involved 
in the reactions. The energies of the colliding hadrons are assumed to be very 
high while the transverse momentum of the partons is low. They all propagate 
in the same direction. When calculating cross sections the colliding particles are 
treated as free particles. This approximation is only valid since the interaction 
between partons is weak for high momentum transfers, i.e. the coupling as{Q'^) 
decreases with increasing momentum transfer Q^. This behavior of a<j(Q^) is 
called asymptotic freedom. A short discussion on that can be found in 

If calculations are performed in the infinite-momentum-frame partons carry 
the momentum fraction x of the hadronic momentum. In this frame transverse 
momenta and masses of the partons can be neglected, since they are small com- 
pared to the partonic momenta. Therefore from momentum conservation one 
obtains 

i i 

for both hadrons, where A denotes the hadron including parton a^. Then for the 
kinematics on parton level one finds that 

S = {Pa+ Phf = {XaPA + XbPB? ~ 2XaXb{PAPB) , (2.17) 

where the parton masses were neglected. 
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Figure 2.1: The parton distribution functions of the proton as measured by the 
CTEQ coUaboration in LO (CTEQ6L1) 0]. 



For calculations of hadronic cross sections it is crucial to know the distribu- 
tion of partons in the colliding hadrons. Since there is no way to find them by 
theoretical considerations or perturbative calculations, due to confinement of the 
quarks and gluons, the Parton Distribution Functions (PDFs) have to be mea- 
sured. In this thesis we use the parton distribution functions for the proton from 
the CTEQ collaboration shown in fig. l2.1[ One finds that u- and d-quarks have 
a very characteristic shape with a larger amount of higher momentum fractions 
X. The distributions of all other flavors and the antiquarks have a shape more 
similar to the gluon PDF, since they arise from next-to-leading-order effects in 
loops attached to the gluons as so-called sea-quarks. 

Hadronic cross sections are finally obtained by convoluting the partonic cross- 
sections with the PDFs: 

1 1 

a{S,Q^)= [dxijdx2f^'Hxi,Q^)f-Hx2,Q^)a,j{xiX2S,Q^). (2.18) 

^1 ^2 

Here the partons i and j carry momentum fraction xi and X2 of the incoming 
hadrons hi and /i2- The sum is taken over all parton flavors, except the top quark 
which is not included in the PDFs of the proton. Q denotes the factorization scale, 
usually chosen equal to the renormalization scale. We choose Q to be the average 
mass of the final state particles. While we are able to do the partonic integration 
in an analytic way, all integrations on the hadronic level obviously have to be 
performed numerically. Therefore we make use of the Monte-Carlo integration 
routine Vegas [lO]. More details about the parton model can be found in [11|. 
Details on our hadronic calculations can be found in the chapters H] and O 

To calculate the PDFs at any given scale, the running of ag has to be consid- 
ered. To be consistent with our choice of the leading-order PDFs, we consider the 
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one loop running of Ug with five quarks, and a given value of 

Ag),j^ = 0.165 GeV . (2.19) 

The dependence of the renormalized strong coupling constant as(/i) on the renor- 
malization parameter /i in the MS-scheme is described by the renormalization 
group equation 

din ij? 

with 



fi{as) , (2.20) 



and 



/?(«,) = -^/3o (2.21) 
47r 



/3o = 11 - ^Nf . (2.22) 



Nf denotes the effective number of quark flavors, i.e. the number of quarks with 
a mass much smaller than the energy scale // of the process. The solution of 
Eq. ()2.20p at one-loop order reads 

a,{^) = -M^o) , (2.23) 

l + ^a.(^o)lng' 



often also noted as 



a. 



4^ /3oln^ ' 

"7 



(2.24) 



with Anj: accordingly defined. If one wants to match the top-quark and other very 
massive particles from SUSY to the running of a<j, one usually requires 

af\mtop) = af\mtop) (2.25) 

which ensures steadiness of the running a<j. At a scale of 580 GeV, which is the 
average mass of the final state particles for our process under investigation, and 

(5) 

including five quark flavors and the given Aq^j-,, one flnds Og = 0.100375. Other 
heavier particles from supersymmetry could, in principle, also appear in loops. 
Due to their higher mass their effect on should be small at LHC energies. If 
SUSY is studied at energies much higher than the masses of all supersymmetric 
particles, the whole particle content of SUSY should be included in the running, 



as it is explained in 1^ . In this case one flnds Pq = 3 for the running of the QCD 
coupling. 
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Until today the Standard Model of Elementary Particle Physics (SM) is very 
successfully in describing the observed phenomena in particle physics. But from 
theoretical considerations we expect that the SM is not the fundamental theory 
of particle physics. 

The SM describes the electroweak and strong interactions of elementary par- 
ticles based on the gauge group SU{3)c x SU{2)l x U{1)y- The gauge group 
SU (3)c describes the strong interaction between quarks and gluons. This theory 
is called Quantum Chromodynamics (QCD). Electroweak interacting fields are 
described by the SU{2)i x U{1)y gauge group. The SM contains 21 parameters 
which can not be derived from fundamental principles and have to be measured. 



Reviews of the Standard Model can be found in e.g. [a], [13[]. An obvious weakness 
of the SM is that gravity is not included as one of the fundamental forces. String 
theories seem to be the best candidates for this. 

Though the SM agrees with the observed phenomena and all predicted par- 
ticles, except for Higgs boson, are found, it can only be an effective low-energy- 
theory since at the Planck scale Mp ~ 10^^ GeV quantum gravitational effects 
play an important role. The question why the ratio ~ 10^'' is so large can not 
be answered by the Standard Model itself and appears quite unnatural. This is 
called the hierachy problem. It results in the instability of the energy scale of the 
Higgs boson when higher order corrections of the renormalizable Standard Model 
are calculated. The mass of the Higgs boson 

mn = v]l^ (3.1) 

depends on the vacuum expectation value v = 246 GeV and the strength of the 
Higgs self-coupling A from the Higgs potential 

y = V0V + ^(</'W, (3.2) 

with A > and /i^ > 0, which is essential for the spontaneous symmetry breaking 
mechanism. At one loop order, the 4-boson interaction in eq. (j3.2|) yields 

AA^^V (3.3) 

in addition to the —fi'^(p^(p term in the potential V, where A is the cut-off param- 
eter. The parameter A represents the energy-scale where new physics appear and 
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the Standard Model looses its validity. These corrections to the mass appear in all 
renormalizable theories, i.e. a theory without divergences after renormalization. 
They are worst for scalar particles. Here one has to choose the parameters A, A 
and consistent with the phenomenologically fixed masses and the vacuum expec- 
tation value. This is usually called fine tuning. A more natural way to stabilize 
the scale of the Higgs mass would be to find a symmetry that makes the correc- 
tions for scalar particles disappear a priori. This is achieved by Supersymmetry 
(SUSY) which provides a fermionic partner to every boson and the other way 
round. Therefore fermions and bosons propagate in each loop and contribute to 
the mass corrections with different signs, due to the relative minus sign for closed 
fermion loops. Terms which are quadratic in A cancel if the coupling of fermions 
to the Higgs boson is equal to the Higgs self-coupling A. The corrections of the 
Higgs self-energy diagrams, including a fermion and the Higgs itself, yield a term 



which can be of the order of the if the fermion masses are not too large. This 
means that SUSY solves the hierachy problem and stabilizes the scale mn such 
that the corrections do not push it towards the Planck Scale Mp. To achieve the 
desired effect, the new supersymmetric partners can not be much heavier than 
1 — 10 TeV [ij]. This makes SUSY most likely observable at the LHC and leaves 
less room for SUSY if there are no new particles found. 

Today physicists are searching for a theory including all forces and it is gener- 
ally expected that at a very high scale electromagnetic, weak, strong and in the 
end also gravitational forces have the same value for the coupling constant. All 
particles are then arranged in a large multiplet, described by a larger symmetry 
group, interacting by only one force. This can be explicitly calculated by using 
renormalization group equations. While unification of the coupling constants does 
not work for the SM, the agreement with the idea of unification is much better 
for SUSY, due to supersymmetric loop contributions. 

A further problem of the SM is that it does not solve the cosmological dark 
matter problem. Today it is known from various cosmological and astronomical 
observations, e.g. investigations of rotation curves of galaxies, that there has to 
be a large amount of dark matter. A candidate for this dark matter could be 
the lightest supersymmetric particle (LSP) if R parity, a newly introduced SUSY- 
quantum-number, is conserved. As a consequence the LSP does not decay to 
lighter SM particles. 

But SUSY is not the only candidate for physics beyond the Standard Model. 
The dark matter problem could also be solved by Kaluza-Klein-number conserving 
extra dimensional theories. There one assumes that the particles known from the 
Standard Model can also propagate into a fifth or even more dimensions, inspired 
by higher dimensional string theory. In this thesis we will always use the simplest of 
all possible extensions to the Standard Model. In the case of Supersymmetry this 
is the Minimal Supersymmetric Standard Model (MSSM). As an extra dimensional 
model, we refer to a model with only one universal extra dimension. 




(3.4) 
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3.1 Supersymmetry 

In this section we want to describe how a super symmetric Lagrangian is con- 
structed. In SUSY every particle from the SM receives equally many supersym- 
metric partner-degrees-of-freedom with the same quantum numbers except for the 
spin. The generators of a supersymmetric theory must turn fermionic states into 
bosonic states and vice versa: 

Qlboson) = |fermion) , Q|fermion) = |boson) . (3.5) 

In general there can be more than = 1 independent supersymmetric operator 
but for phenomenological studies we will assume the easiest case of = 1. If SUSY 
is an exact symmetry, the masses of particles and their supersymmetric partners 
would be equal. But since the SUSY partners of the known SM particles are not 
yet found, SUSY has to be broken. 

The generators of the SUSY transformations, as given in eq. ()3.5p . satisfy the 
following formulae: 

{Qa,Qb} = 0, 

{QlQl} = 0, (3.6) 
{Qa,Ql} = a'^.P,. 

At this point one can already see that SUSY is even more than just a symmetry 
between fermionic and bosonic degrees of freedom but also an extension of the 
usual space time symmetry, since the momentum four vector is linked to the spin 
of the particles. 

In SUSY all particles and their partners are arranged in two different kinds 
of supermultiplets. One is the matter or chiral multiplet that includes a two- 
component Weyl fermion and two real scalars called sfermions. The other is the 
gauge or vector multiplet including a spin 1 vector boson and a spin ^ Weyl 
fermion called gaugino. Formulating a supersymmetric theory is possible by using 
the superfield notation or the more familiar spinor notation. 



Following the discussion of iJ], we introduce the easiest "supersymmetric" 
Lagrangian in the two-component spinor notation, including only one complex 
scalar spin-0 field (f) as well as one L-type, i.e. a left-chiral spinor field x- 

c = d^<p^d^'<p + x^iaf'd^x ■ (3.7) 
The fields <j) and x are linked by the following SUSY transformations: 

Sa = -[ia^{ia2n]M ■ (3-8) 

Here ^ is not a field, because it is independent of x, but a constant L-type spinor 
parametrizing the SUSY transformation. The transformations for their hermitian 
conjugates are found to be 

5^^t = x\i<T2)C, 

= -d,,<P^^^ia2iaf' . (3.9) 
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Under these transformations the Lagrangian £ changes only by a total derivative, 
leaving the action and the equations of motion invariant. But in order to show 
that the theory is supersymmetric, one also has to show that the SUSY algebra 
closes, i.e. that a term like {Se^^ — S^St)4' with (p being a fermionic or bosonic 
field of our chiral multiplet vanishes. Unfortunately this does not hold for off-shell 
particles since the Weyl-equation does hold and therefore does not take out two 
fermionic degrees of freedom of the L-type spinor. As a result one has two bosonic 
and four fermionic degrees of freedom in the complex fields (p and x- Due to the 
fact that in every multiplet the number of degrees of freedom has to be equal, 
we have to introduce two more by a second scalar field F. The Lagrangian then 
becomes 

/:chirai = d^^^d^'^ + x^ia^'dx + F^F . (3.10) 
The SUSY transformations leaving this Lagrangian invariant are: 

6^F = -i^^a^'d^x, 
5^F^ = id^x^a''^, 

= -[icT'^(icT2r)]9;.</' + Ft^^ (3.11) 

SaX^ = -d^,ct)^(^ia2ia>' + iF . 

Under these transformations the supersymmetry algebra is closed. The first su- 
persymmetric theory with interacting particles was formulated by Wess and Zu- 
mino [H. They added the most general and renormalizable term 

Ant = '^^)^^ - • ^'J + ^-C- • (3-12) 

Certain conditions for the superpotential W can be derived and after eliminating 
the non-physical degrees of freedom, i.e. the scalar field F, one finds 

C-^NZ = Uee " \Wi\^ - ^{Wij Xi ' Xj + h.C.} , (3.13) 

with 

W = ^Mij4>i<f>j + ^yijk4>i4>j(t>k , (3.14) 
dW 

Wi = j^, (3.15) 

= ——. (3.16) 

W is called superpotential. Only taking into account terms which are quadratic 
in the fields and calculating their equations of motion one finds that for the fields 
(/> and X the usual Klein-Gordon-equation is satisfied for each component. Both 
fields have the same mass M coming from the interaction term in the Wess- 
Zumino model. This shows that in unbroken SUSY masses are equal within the 
same supermultiplet. 

Beside the chiral multiplets one also has to take into account the vector or 
gauge multiplets of the different gauge groups. Here we meet the same problems 
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concerning tlie degrees of freedom as before and we liave to introduce a second 
auxiliary field. For the SU{2) and SU{3) gauge groups the only difference is 
located in the different generators of the gauge groups and therefore in the number 
of corresponding gaugino fields. The resulting supersymmetric Lagrangian for the 
gauge multiplet is given by 

>Cgauge = -If^.F'^'"^ + iX'^^a'^iD^Xr + ^D'^D'' . (3.17) 

Here is a scalar auxiliary field, contributing the additional degrees of freedom, 
F^"^ is the usual Maxwell field strength tensor and A is a L-type spinor, the gaugino. 
In the case of an SU(2) gauge group a triplet of gauginos is needed while for the 
SU(3) gauge group a gaugino octet appears. The same applies to the auxiliary 
field D". For the SUSY transformations one finds 

S^D" = -i{^^a^{D^Xr-{D^Xpa^^O- (3-18) 

In order to retain gauge invariance of the supersymmetric Lagrangian one has to 
replace the usual derivative by the covariant derivative. 

Since we now have a SUSY invariant gauge and chiral supermultiplet La- 
grangian we can ask for interactions between particles from those two multi- 
plets. After dimensional considerations, inserting all possible interaction terms 
and checking SUSY invariance and gauge invariance, one obtains the following 
Lagrangian for the combined gauge and chiral multiplets in the non-abelian case: 

■^gauge & chiral ~ -^gauge "I" -^chiral "I" -^int. gauge & chiral 

= - ^Fj^^F'"'"^ + iX'^^ai'iD^X)'^ + ^D'^D'^+ 

+ D^4D^'<Pi + xl^^^D^Xi + F^Fi+ (3.19) 

_d(j)i 2d(j)i(t)j 
-V2g [i^T^Xi) ■ A" + A-t • (xfr^,)] - 5(4^^"</'^)^" ■ 

Both parameters ^ parametrizing the SUSY transformations of chiral and gauge 
multiplet do not necessarily agree. In fact one finds for the gauge multiplet that 
the parameter ^ has to be changed into — ^ / \/2. The transformation of F- for the 
combined Lagrangian has to be changed and is now given by 

6^fI = -V2g4T''C ■ A" + ^^^x^a^^C . (3.20) 

Since the fields F and D can be replaced in the given Lagrangian, one obtains the 
whole scalar potential consisting of the so-called F- and D-term, as 

y(0,,4) = + ^EEE^G {4ts4>-^) (^jrg^,) , (3.21) 

G a i,j 
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where in general the sum is taken over different gauge groups, their generators 
and the arising fields. 

There are even some more possible interaction terms which keep the La- 
grangian invariant under SUSY transformations and are renormalizable. But they 
either violate lepton number L or baryon number B. Since one did not observe L 
and B violating processes, e.g. proton decay, one has to exclude these terms from 
the Lagrangian by claiming a symmetry to hold. This is the so-called i2-parity, 
given by 

R = {-lfB+^+^\ (3.22) 

where 6 is 1/3 and —1/3 for quarks and antiquarks, L equals 1 and —1 for leptons 
and antileptons and s corresponds to the spin of the particle. Therefore R equals 
+1 for conventional particles and —1 for their SUSY-partners. As a consequence 
SUSY particles can only appear in pairs at a given vertex. In our calculations we 
assume i?-parity to be always conserved. So the lightest supersymmetric particle 
(LSP) is stable and therefore a serious candidate for dark matter. This LSP is 
usually the lightest neutralino. When SUSY searches are performed at the LHC, 
two times the mass of the LSP is the minimum amount of missing energy. 

3.1.1 The MSSM 

Since we now have an understanding of how the multiplets are set up, we present 
the particle spectrum of the MSSM in table [3Tl Due to its length we do not want 
to introduce the Lagrangian for the physical fields. It can be found in [IB] and [l^ . 
In the last section we mentioned the superpotential W which has not been specified 
yet. The MSSM is given by the following choice for the superpotential 

W = vi^UiQjH^ - vi^diQjHd - yl'e^LjHd + fiH^H^ , (3.23) 

where the contraction over gauge indices is implicit. A new free parameter fi, the 
Higgs self-coupling is introduced in this superpotential. Q, L, e, u and d denote 
the superfields as they are given in table 13. 1[ and have to be understood 
as 

For anomaly cancellation it is important to have two Higgs fields, both having 
opposite hypercharge. In order to break SUSY one introduces so-called soft SUSY- 
hreaking terms and breaks the symmetry explicitly. These terms are called "soft" 
because they do not destroy the solution of the hierachy problem by reintroducing 
new divergences. 

The origin of the soft SUSY-breaking terms is far from clear. Explaining the 
different possibilities for their existence is non-trivial and goes far beyond a short 
introduction. Including the SUSY breaking terms more than 100 parameters are 
added to SUSY, making concrete predictions quite difficult. When a special SUSY 
breaking mechanism is assumed this number is drastically reduced. The mecha- 
nisms most often discussed are minimal supergravity (mSUGRA) , gauge- mediated 
SUSY breaking (GMSB) and anomaly-mediated SUSY breaking (AMSB). SUSY 
is broken in a hidden sector in all the scenarios and is mediated to the visible 
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Table 3.1: Particle content of the MSSM with three families for quarks and leptons. 



sector by gravitational or gauge interactions or by the Super- Weyl anomaly. A 



short overview of these SUSY breaking schemes is given in [17|] . 

In the mSUGRA scenario the whole parameter space is described by four 
parameters and a sign. These are the scalar mass parameter ttiq, the gaugino 
mass parameter mi/2) the trilinear coupling ; the ratio of the Higgs vacuum 
expectation values tan/3 and the sign of the Higgs sector parameter \i. From these 
parameters one obtains the mass spectra of the SUSY particles by renormalization 
group running from the SUSY breaking high-energy scale to the weak scale where 
they can be used for calculations of cross sections or decay branching ratios. 

But even for a number of five parameters, instead of more than 100, scanning 
the whole parameter space for simulations of SUSY signatures is much too time- 
consuming. Therefore one uses selected parameter points [l^, which are consistent 
with all constraints from phenomenological studies, like the dark matter content of 
the universe and SUSY searches at Tevatron. For example, a region with tan/3 < 3 
is excluded by LEP experiments. 

As our standard scenario we choose the parameter point SPS la, which is a 
typical mSUGRA scenario. It is specified by 

mo = lOOGeV, mi/2 = 250GeV, = -lOOGeV, tan/3 = 10, /i > . 



3.1.2 How to treat Majorana Particles 

Unlike in the SM, in SUSY various Majorana particles appear. Majorana particles 
have the interesting feature that they decay equally often into particles and an- 
tiparticles. Therefore processes like the gluino pair production can be used to find 
the helpful like-sign dilepton signature in order to claim that SUSY is observed. 
This is especially interesting at hadron colliders, since the production cross section 
for gluino pairs is quite large [§]. 
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Compared to off-shell fermions, the Dirac equation reduces the number of 
degrees of freedom of on-shell fermions by a factor of two. Off-shell Majorana 
fermions have, instead of eight degrees of freedom as usual fermions, a reduced 
number of four degrees of freedom. This is obvious from the general condition for 
Majorana fermions 



Therefore they can be represented in a unique way by two component complex 
spinors or x as well as by a four dimensional spinor 

K = ( J^^^. ) or equivalently = ) . (3.26) 

The spinors ip ox x are defined by using the chirality projection operators as 

Pr^= (q) and Pl^= , (3.27) 



Pr=CA?-) and Pl={^-^]. (3.28) 



with 



In [_19|] Feynman rules for Majorana fermions are derived and it is proved that 
it is possible to treat them nearly equal to Dirac fermions. Therefore one chooses 
an arbitrary direction for the fermion flow and writes down the Feynman rules 
for the Majorana fermions in the same way as one does for Dirac fermions, going 
against the fixed flux. For the external Majorana fermions one uses the same 
spinors u(p, s), i;(p, s) and the same spin sums 



u'^ (p)u'^ (p) = p + m and ''^^ v'^ (p)v'^ (p) = p — m (3.29) 



as in the case of Dirac fermions. For the correct choice of the external spinor, 
depending on momentum and direction of the fermion flow, we use the Feynman 
rules given in [l9|]. 

Of course Majorana fermions have no arrow on their lines indicating if they 
are particle or anti-particle since they are both. These arrows, belonging to Dirac 
fermions, indicate the so-called fermion number flow. 

Since one fixes the fermion flow in the Feynman diagrams there are two dif- 
ferent Feynman rules for a vertex containing Dirac and Majorana fermions, one 
with the fermion number flow parallel and one anti-parallel to the direction of 
the fermion flow. These two rules are, as in the case of the gluino-quark-squark 
vertex, associated with the helicity of the squark: 

Cqgq = -V2gst% ^ {la Pl q} q*Lf + q) Pr ga q{f - la Pr q) qtf " <ff Pl 9a 4/) • 

f=u,d 

Here g denotes the Majorana gluino, q a usual Dirac quark and q a scalar squark. 
Obviously so-called left-/right-handed squarks only couple to left-/right-handed 
quarks. The Feynman rules for SUSY-QCD, including quark-gluino-squark inter- 
actions, can be found in the appendix. 
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3.2 Extra Dimensions 



Though Supersymmetry is the most frequently studied extension of the SM, there 
are also other theories solving problems like the dark matter and the hierachy 
problem. One of them is similar to the Standard Model extended by additional 
space-like dimensions, called Extra Dimensions. There are different kinds of extra 
dimensional theories discussed in the literature, beginning in 1919 when Kaluza 
and Klein unified gravity and electromagnetism using a hidden fifth dimension, 
compactified on a circle [i^, [2]]]. Later it was realized that string theories can 
not be formulated in four dimensional space-time. In the following years, different 
ideas like flat, large extra dimensions on a millimeter-scale were proposed in |22!] in 
1998 (ADD). In this model only gravity, i.e. the graviton, can propagate into the 
extra dimension, while all SM particles are confined to the four dimensional space- 
time. In contrast to this, a theory was proposed in 2001, where all SM particles 
can propagate into the extra dimension [2] . This kind of extra dimensional theory 
is therefore called Universal Extra Dimension (UED). In the literature one can 
additionally find calculations with fermions confined to the usual 4D world, while 
gauge bosons can also propagate into the extra dimension 23|], 24]. 

Various problems can be solved by UED models. Including boundary condi- 
tions the symmetry breaking in the electroweak sector can be performed higgsless, 
being an interesting option if Higgs searches at the LHC fail [is^] . There is also 
a solution for the dark matter problem of cosmology. When Kaluza-Klein (KK) 
parity, a quantum number being quite similar to the R-parity in SUSY, is con- 
served, a lightest electrically uncharged and colorless KK-particle is present. This 
particle is an interesting dark matter candidate [2^. 

It is known that higher dimensional operators violate global symmetries like 
baryon number and lepton number and can create flavor changing neutral currents. 
The SM, by chance, has no problems with baryon-number violating interaction 
terms mediating the proton decay, since all possible terms are of higher dimension 
and therefore non-renormalizable. In contrast UED in principle can include these 
higher dimensional operators. But since the proton decay is limited to very low 
decay rates by experiment, there should be a mechanism suppressin g th ese higher 
dimensional terms in UED. Therefore propositions like split fermions [27[ or simply 
a suppression by the cut-off scale are made. 

In a model with universal extra dimensions, i.e. all particles of the theory 
can propagate into the extra dimension, the size of the extra dimension is of the 
order of ii ~ TeV~^. Otherwise their existence would contradict observations 
since they have not yet been found. This can only be solved by compactifying the 
extra dimension on a scale where large energies are needed to see their effects. In 
our calculations we use this UED model and only one additional space-like extra 
dimension. 

Here we concentrate on the collider implications of UED. We do not want to 
give an overview of the different extra dimensional theories but instead explain 
in detail how to derive the Feynman rules and the UED Lagrangian in four di- 
mensions. A collection of UED Feynman rules can be found in the appendix. For 
simplicity we concentrate on the QCD part of the Lagrangian and derive the cou- 
plings of the neutral KK-gauge bosons. In this thesis we do not intend to show that 



20 



Chapter 3 Standard Model Extensions 



UED is the theory reahzed by nature but contrast it to SUSY to investigate the 
differences coming from the different spin assignments. While the SUSY particles 
have different spin compared to their SM partners, particles in UED have the same 
spin statistic behavior as their SM partners. For example, gluons have bosonic 
partners in UED (heavy KK-gluons) while they have fermionic partners in SUSY 
(gluinos). One can find numerous papers on how UED could mimic the signature 
of a supersymmetric theory, e.g. [28|], though SUSY and UED have fundamental 
differences. While usually the mass spectrum in UED is quite degenerate and 
thereby different from SUSY, radiative corrections and boundary conditions can 
make the spectrum more SUSY-like. Reviews on extra dimensional theories can 



be found in 271] , [29| and in [30|] . The latter one focuses on collider phenomenology 



and dark matter. 

3.2.1 Deriving a UED-QCD Lagrangian in 4D 

In this chapter we derive the 4D effective Lagrangian and the UED-QCD Feynman 
rules. The basic idea is to generalize the 4D Lagrangian to its analog in 5D, then 
to Fourier expand the Lagrangian in the coordinate y of the extra dimension. 
Orbifold compactification of the fifth dimension is performed and the usual SM 
terms plus additional interactions of KK-excitations are obtained, i.e. one obtains 
the 4D Lagrangian by integrating out the extra dimension. From the Lagrangian 
we get the masses and couplings of the new particles. Since our discussion is 
parallel to those in 3, and H, 

we denote the multiplets of the Standard 



Model by 

Qf*^(x), Ufi^Hx) and Df^{x) . (3.30) 

As we will see later, each of the SM particles has infinitely many KK-partners. 
These are equidistant in mass and form a KK-tower. The different excitation 
levels of a tower are denoted by an index n. At the LHC, we can only expect the 
first KK-partner of this tower to be produced if R is sufficiently small. Therefore 
in our calculations we do not take into account the heavier KK-partners. 

Since SM quarks are assumed massless in our calculations they can be denoted 
as two component Weyl-Spinors in 4D. In the 5D theory quark multiplets consist 
of massless four component quark fields Q{x, y), U{x, y) and D{x, y). There are no 
chiral fermions in a 5D theory, since it is not possible to construct a 7-matrix in an 
even number of spatial dimensions, having the same anti-commutativity property 
as 75 in 4D. Bilinears containing 75-matrices are not Lorentz-invariant in a 5D 
theory and can therefore not be included in the Lagrangian. 

The additionally needed 7-matrix is chosen to be 

^ (_^)D/2+s+l ^0^1 _ _ ^ (3 

where 7^^, . . . ,7^"^ are the 7-matrices in even dimensions D and s is the number 
of space-like dimensions. Since this 7-matrix anticommutes with 7^^, . . . ,7^"^, it 
can be used as 7^ in five dimensions. This choice for the additional 7-matrix is the 



usual choice in a higher dimensional supergravity theory, as explained e.g. in 31 1 
In the case of a 5D theory this yields 
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For 7^ = the Chfford Algebra 



lixlu + Ivlti = 2g^^ = 2 diag(+l, -1,-1,-1,-1) 



(3.33) 



still holds but it does not anticommute with 75 as the other 7-matrices 7°, . . . , 7^. 

We are forced to introduce two 5D fermionic fields with the quantum numbers 
of a left- and a right-chiral spinor. Decomposing these 5D fields into 4D one obtains 
a left- and a right-handed zero-mode. Of course, since right-handed doublets and 
left-handed singlets are not observed, at least their zero-modes have to vanish. 
Therefore one chooses a S1/Z2 compactification. Higher modes n ^ 1 can, in 
general, contain both, left- and right-chiral particles without a contradiction to 
experiments. In order to identify Qf^\x),U^^ {x) and D^\x) with the SM 
particles, the zero-modes of the Fourier expansion of Q{x,y), U{x,y) and D{x,y) 
must be even under y ^ —y. This yields the following Fourier expansions for the 
doublet Q and the singlets U and D: 



Q{x,y) = 
Uix,y) = 
D{x,y) = 



1 
1 



u{x) 
d{x) 



+ 



00 

V2 E [qUx) cos (^) + Qlix) sin (^) 



n=l 



00 

ur{x) + V2 E [U^ix) cos (^) + U2{x) sin 

n=l 

00 

dnix) + V2Y^ [dUx] cos {^) + Dlix) sin 



n=l 



Since the zero-mode of a gauge field, polarized along a direction in 4D, must be 
even under y ^ — y it must be odd when it is polarized along the y axis. Otherwise 
we would find unphysical massless zero-modes in the 4D effective Lagrangian. 
Therefore its Fourier expansion is given by 



1 



A^^,{x) + V2Y,A%n{x)cos {^) 

n=l 

V2 °" 



At{x,y) 



I-kR 



n=l 



(ny\ 
x)sm(-j . 



(3.34) 
(3.35) 



In unitary gauge the last term disappears, due to the gauge choice A\^(x) = 0. 
Normalization of the zero modes is different from the higher modes in order to 
obtain canonically normalized terms in the kinetic part of the Lagrangian and 
because of the limits for the integration of y. The Lagrangian in 5D is given by 



£5 = iQ{x, y) {d^ + zffst" Al,{x. y)\ Q(x, y) 



(3.36) 



with g = ;^^=, where g equals the usual SM coupling constant gg and R is the 
radius of compactification. The index m is equivalent to the usual 4D index /x, 
here also including the extra dimension, i.e. m G {/x, 4}. The KK-mass eigenstates 
are given by 

QX^{X) = Pl,R ^^'ll^j ) , \J\d^) ^ Pr,l UnM^) , D^j^ix) = Pr,l dn,2{x) , 
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with Pifi = ^(1 =F75)- The fields and qn^2 denote two towers of KK-partners, 
arising for each usual SM quark q. Here the first index denotes the excitation 
level of the KK-tower while the second index denotes the tower itself. Each of the 
fields Q,U and D gets a contribution to its mass from integrating the kinetic term. 
This is different from the SM, where fermionic masses only come from the Yukawa 
interaction terms and the electroweak symmetry breaking. In principle the fields 
Qn^i and qn,2 could also get contributions to their mass from symmetry breaking. 
Integration of the kinetic terms, after insertion of the Fourier expansions of Q,U 
and -D, yields the effective 4D kinetic and mass terms 

ttR „ , 

u{x) 



ij Q{x,y) 



■y"^dmQix,y)dy = i 



[u{x)d{x))Li^d^ 



d{x) 



n=l 

71 — 71 — 

+ i^Ql{x)Ql{x) + i-Ql{x)Ql{x) 



(3.37) 



i I U{x, y) j'^dmUix, y)dy = i [(n(x))^7^9^ (n(x))^ 



oo 

+ Y.^R{^)^'W{^) + ui{x)rd,ui{x) 

n=l 

71 — 71 — 

+ i^Ul{x)Ul{x)+t-Ul{x)Ul{x) 



(3.38) 



and similarly for the second singlet D{x,y). The last term, i.e. the mass term 
mixing left- and right-handed fields, derives from integrating the extra dimensional 
term with m = 4. The 74-matrix disappears from the Lagrangian since 



75(1 =F 75) ==F(1=F75) • 



(3.39) 



One also receives the mass of the heavy KK-partners M„ = ^, depending on the 
size of the extra dimension. Contributions to the KK-masses from the Higgs part 
of the Lagrangian are neglected here. The interaction term is the second term in 
our Lagrangian in eq. (j3.36p . given by the usual covariant derivative, and yields 

nR 

-9,j Q{x,yh"'eA-^{x,y)Q{x,y)dy 




g-L(x)7^t"?L(x)A^,o(^) + E [Qli^h't^'Qli^) + Qirt'^Qlix)] A;^oi^) 

n=l 

00 

+ E i^Lixh^f^Qlix) + Qlixh'-eqLix)] A^,„(x) (3.40) 

n=l 
^ 00 

+ 71 E [Qli^h^i''QT(.^m,\n.-n\+Si,„,+n) 



n,m,i=l 

Ql{x)^n-Q^{x){5,^\^_^\ - 5,,„+„)] A^i 
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For U and D fields we again obtain a similar expression given by 
-9, J U{x,yh^eA''^{x,y)U{x,y)dy 





n=l 

oo 

+ [uR{x)rt''U^ix) + U^{xhH-UR{x)] (3.41) 

n=l 
^ oo 

n,m,£=l 

+ U2{xh^t'^Ur{x){6e,lm-n\-Sl,„,+n)] 

Expressing everything in terms of the mass eigenstates g^.i and qn,2, and summing 
up the contributions of both heavy quark towers, one obtains 



:,m-\-n) 



int 



-g\q{x)^H-q{x)A1^^{x) 

oo 

+ Y ['ln,l{xh''t''qn,l{x) + qn,2{^h^t^(ln,2ix)] A^^^oi^) 
n=l 
oo 

+ Yl ['lL{x)j^'eqn,l{x) + qn,l{xh^t\L{x)] Al^^{x) 

n=l 
oo 

qn,2{x) + qn,2ix)j^t°'qR{x)] j4"„(x) 

n=l 
^ oo 

+ ^ [-^n,i(a;)7''75i"gm,i(x) + g„,2(x)7''75t''gm,2(a;)]^;^,^5£. 



V2 
1 



(3.42) 



m+n 



n,m,l=l 
oo 



l,\m—n\ 



+ 7f 2^ [g„,i(x)7'^t'^g™,i(x) +g„,2(x)7''t"gm,2(x)]^;i,^<5. 

n,m,£=l 

Due to the field expansions, integrations include sine and cosine integrals, resulting 
in the Kronecker-(5s. Interactions between different towers, qn^i and qn,2^ which 
are a mixture of the singlet and doublet particles, are automatically excluded for 
symmetry reasons and KK-number is automatically conserved. The relation 

I ± ± • • • ± nj^_^ \ = , (3.43) 

deriving from the Kronecker-(5s, implies KK-number conservation. But it is only 
valid at tree and broken at loop level. 

To show in an example that KK-number is conserved at each vertex we consider 
the vertex of two left handed fermions coupling to one gauge boson. This is one 
part of the Lagrangian already derived in eq. (j3.40p . Therefore one has to integrate 
out the 5D Lagrangian 



9, / Ql{x, y)7"^t"A^(x, y)QL{x, y)dy . 



(3.44) 
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Inserting the Fourier expansion of the fields, the term including only the zero- 
modes of the fermions reads 



tR 



-9, 



1 







n=l 



QLix,y)dy . 



Since integration of the cosine gives 



I 





cos ] dy = ttR 5n,o 



(3.45) 



the zero-modes of the fermions only interact with the zero-modes of the gauge 
bosons. Considering the interaction between one zero-mode fermion, excitation 
level / of a second fermion and the gauge boson, one finds 



9. 



qL{x,y)^"'e- 



1 



1=1 



(3.46) 



X 



^V2^QUx)cos[^)dy 

n=l 



Now integration of the cosines yields 

ttR 

cos (- j cos 



dy = vri? 5. 



n,l 



(3.47) 



Here it becomes obvious that no mixing between different KK-excitation levels 
occurs since there is only interaction if n = /. The last interaction term derives 
from two KK-fermions and one KK-gauge boson. Therefore we integrate eq. p.44p 
after inserting all Fourier expansions and find the integral 



ttR 



cos 



^) ("^) (^) '^^ ^ ^ "^'.l™-"! + ^l,rn+n , (3.48) 



where I is the KK-level index of the gauge boson and n and m are the indices 
of the excitations of the fermions. We again find that KK-number is conserved 
also for this vertex. Of course, we can easily violate KK-number at loop level by 
choosing one vertex with a zero-mode scattering into two first level excitations 
and the other vertex with the two first level excitations creating one second level 
excitation. But nevertheless KK-parity, given by 



pkk = i-iy 



(3.49) 



is still conserved at the loop level. If one assumes that the SM fermions can not 
propagate into the extra dimension and therefore multiplies the Lagrangian in 
eq. (I3.44P by d{y), the SM fermions couple to all gauge boson excitations with 
the same strength. Since the absence of the (5-distribution results in KK-parity 
conservation, this can be understood as a consequence of translation invariance in 
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the fifth dimension. This conservation is the reason for the lowest lying KK-modes 
to be stable and thereby an interesting candidate for dark matter. It also implies 
that at colliders, KK-particles can only be produced in pairs. 

We also want to derive the purely gluonic part of UED-QCD which we have 
not taken into account yet. In equivalence to the Standard Model, the purely 
gluonic part is given by 



' ^ ^ mn ^ 



with the field strength tensor in 5D given by 

■pa _ p. Aa _ ^ Aa _ rabc Ab A' 
mn ~ '-'m.-'^n '-'n-^-m HbJ 



(3.50) 



(3.51) 



and m,n £ {fx, 4} are the 5D space time indices. The first term in brackets in 
eq. (|3.50p contains interactions and kinetic terms of gluons g and their heavy 
partners 5*, given by 



4 /^^ 



dyAld^'A'"' 



d^Ald'^A^''' + dyA^d^'A'"' 
^gJ^-^^A^^Ali^d^A"" - d^Af"" 

2 rabc pade j^b j^c 
t/5 J J fj, 1^ 



(3.52) 



while the second term in brackets in eq. (I3.50p becomes the heavy gluons mass 
term 



(3.53) 



where we again apply the gauge A'^ = 0. 

Now integration over the Lagrangian has to be performed after inserting the 
expansion (j3.34p . For the three-gluon interaction term one obtains 



1 



■nR 

;g,r''l Al{x,y)Al{x,y) d^^ A'^'^ {x , y) - d'' A^^^ {x , y) 




dy 



■\9f 



A%{x)Al,{x)[d^A'(,'^{x)-d''A^,%x 

00 

+ ZA%{x) AUx) [aMr(x) - d^A^^ix 

n=l 

+ 71 E A'^MKm[d''Ar{x)-d^A^'{x) 

n,m,£=l 



(3.54) 



'.itmitn 



The four gluon interaction in 4D is given by 
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2 fabc fade 

4 



\gl r'T"" I y)At{x, y)A^''{x, y)dy 

A%{x)Al,{x)A^;'{x)A^,'^ix) 

oo 

+ 6A%{x)Alo{x) (x)^r(^) (3.55) 

n=l 

4 

+ -T^^Ui^) Yl AlJx)A';^{x)Ar{x)Si,±m±n 

n,m,e=l 

^ oo 

+ 2 E A'^^{x)Al^{x)Af{x)Ar{x)6k,±m±n±i . 

n,m,l,k=l 

The Kronecker-5s again derive from integrations over sine and cosine functions 
and have to be understood in the fohowing way: 



•,k+m + Se,k-m + Se,m-k + ^l-m-k ■ (3.56) 



The last Kronecker-(5 can never be non-zero, since n, m and / are defined as positive 
integers. 



The heavy gluons mass term from eq. ()3.53p gives 

- 2 / d,A-ix,y)d'A>^-{x,y)dy = --^ ^ ^«„(x)^r(^) • (3-57) 
"=i 

For the KK-gauge bosons we again find M„ = ^. This leads to the naive con- 
clusion that the mass spectra within one mode in UED are highly degenerate, 
neglecting other symmetry-consistent boundary terms and the Higgs mechanism. 
This means that the mass spectra, which are typically expected in UED, are quite 
different from a typical SUSY spectrum, e.g. the SPS la benchmark point. 

All Feynman-rules and Lagrangians were re-derived for UED-QCD. Feynman- 
rules can be found in the appendix. They agree with the rules given in 231] 
and H- 

As will be explained in section \b.2\ we include a width in the propagator of the 
massive gluons, denoted by , and the massive quarks, denoted by qn,i, qn,2- This 
Breit-Wigner propagator will be used for the decay chains in order to regularize 
divergences at the pole. 

Instead of the usual polarization sum for massless gluons g, for the massive 
external gauge bosons (7* with n > 1, we use a polarization sum given in the 
unitary gauge by 

J2<n{k,^)eUk,^) = {-9,. + ^y^' ■ (3-58) 



This is equal to the polarization sum for usual massive gauge bosons. It contains 
three physical degrees of freedom. 
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3.2.2 Electroweak sector in UED 

In general it has to be taken into account that all particles with equal quantum 
numbers can mix. This is important for the couplings of the neutral heavy gauge 
bosons from the electroweak Lagrangian. In this thesis we neglected the Higgs 
part in the electroweak Lagrangian. One can in principle add it, though this is 



not necessary since also higgsless models in UED exist |25l |. Taking the Higgs 
terms into account, one obtains a second mass term for the quarks. These two 
mass terms mix singlets and doublets of KK-particles with the mass matrix 



After diagonalization one obtains mass eigenstates with 



Mn = ]J + {MsM? . (3.60) 

The mixing of the weak eigenstates (5!^"(x) and [/'"(x), forming the mass eigen- 
states Un^i and Un^2, is then denoted by 

un,i = cos a(") Q'„"(x) + sin [/' "(x) , 

un,2 = sina(")75 Q'/(x) -cosa(")75 [/'"(x) . (3.61) 

The mixing angle a^"^ between singlets and doublets is given, without NLO cor- 
rections, by 

tan2a(") = ^ . (3.62) 
7! 

Here Q'^{x) denotes only the upper component of the weak doublet eigenstate 
and [/' "(x) denotes the ti-type singlets. The mixing structure for the d-type com- 
ponent of Q'/^{x) and singlet D' "(x) is equal. This results in the mass eigenstates 
dn,i and dn,2- Since it is suppressed by the KK-excitation mass, the mixing be- 
tween singlets and doublets will, except for the top quark, be quite small. The 
mass term from the kinematic part of the Lagrangian is not present for the zero- 
modes. That is consistent with the non-mixing singlets and doublets in the SM. 

The treatment of the QCD gauge part was exemplary for the SM gauge group 
SU{2>)c X SU{2)l X U{1). In the same way we calculate the Feynman rules from 
the electroweak sector of UED, using the following Lagrangian 



= {u'ix,y),d'ix,y))^j"'(g^A:, + Y,gYB, 



L 



u'{x,y) 

L ' V^T^^™ ' '"^^ \d'ix,y) 
+ f7'(x,y)7™ {Ys,ugY B^) U'{x,y) 

+ D'{x, y) (y,,^ gy B„,) D'{x, y) . (3.63) 

Here U' and D' denote the eigenstates of the weak interaction. Integration of the 
fifth dimension is performed as usual. But since for our later calculations we only 
need the uncharged gauge boson fields, we only compute the neutral part of the 
effective Lagrangian in 4D. 
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As a short example we want to give the Lagrangian for the interaction of the 
first level KK-partner of the SM photon with a SM n-type quark and a first level 
excitation of the KK-quark. As given in eq. (I3.63p . the interaction reads 

Of )(x,y)yd,t;gy7'"i?mQf + tJ'^''\x,y)Y,^u9Yl'^BmU'{x,y) , (3.64) 

where the first term derives from a doublet and the second terms derives from 
a singlet interaction. After insertion of the Fourier expansion for the KK-quark, 
compactification and insertion of the mixing matrix of the weak eigenstates one 
finds 

u 7^" QY Bi ^ Yd^u Pl cos a^^) qi,i - uj^ gy Bi ^ Yd^jj Pl sin a^^^ qi^2 
+ u^y^'gyBi^ Ys^u Pr sin a^^^ qi,i - uj'' gy Bi ^ Yg^u Pr cos a^^^ gi,2 . 

This corresponds to the Feynman rules given in the appendix. 

In [s^] the mixing of the electroweak gauge bosons is investigated and, includ- 
ing corrections from NLO, found to be very small for the first excitation mode, 
given by the corresponding mass matrix 



(3.65) 



In our calculations we assume, that the first level KK-excitations ^3^1 and B\ 

(n) 

already are the mass eigenstates since the mixing angle 0^ , defined by 

7n = smd^^ Az^n + co's.Q^^ Bn (3.66) 
Zr. = zo^Q^^ A^^n-^^^Q^y} Bn (3.67) 

approaches zero because of the degenerate masses for higher level excitations in 

We explicitly calculated only four vertices for the neutral gauge bosons. A 
complete collection of electroweak UED Feynman rules is given in the appendix 
of 
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Chapter 4 

SUSY- and UED-QCD 
Processes at the LHC 

In this chapter we want to calculate the partonic and hadronic cross sections and 
present the contributing Feynman diagrams for squark and gluino production. 
Since the Large Hadron Collider (LHC) is a proton-proton collider it is necessary 
to calculate the hadronic cross sections, according to the parton model, using the 
PDFs for the inner structure of the proton. Since the presented processes are all 
SUSY-QCD processes, they are very important for hadronic colliders. 

4.1 Partonic Cross Sections 

On the partonic level the following processes contribute to gluino and squark 
production at hadron colliders at leading order: 



qq production 



qq production 
gg production 



qg production 

The incoming momenta are denoted ki and k2 while the outgoing momenta are 
called pi and p2- The given processes exist for both chiralities of the squarks 
denoted by qi and qr which are not given explicitly above. In the case of qq 
and qg production charge conjugated processes also have to be considered. When 
presenting the results for squark-antisquark final states, the processes of quark- 
antiquark and gluon-gluon scattering have to be taken into account. Similarly for 
gluino-gluino production, there are contributions coming from quark-antiquark 



g+9 — ^qi+Qi 

qi+qj — *• qi +qj and c.c. 

qi+qi — ^9+9 

9+9 — '9+9 
qi+9 — ^ qi +9 and c.c. . 
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Figure 4.1: LO contributions to the production of squarks and gluons at hadron 
colUders in SUSY. 

and gluon-gluon initial states. The only initial states contributing to squark- 
squark pair production are two incoming quarks. In the case of squark-gluino 
pairs in the final state, only an incoming quark-gluon pair gives a contribution. In 
our calculations all outgoing squark flavor and chirality states are considered to 
have the same mass while the top-squark is excluded from the final state. Incoming 
top-quarks will always be neglected, since in the PDFs the top-quark content is 
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expected to be approximately zero, due to the high top mass. All the other quark 
flavors, Uf = 5, are treated as massless particles since the energy at the LHC 
{^/s = 14 TeV) is much higher than the mass of the bottom-quark. Therefore for 
on-shell particles one uses kf = and pi = mf . 

Our notation for the matrix elements follows the notation of 32] ■ We use the 
Mandelstam variables which are kinematical invariants and are defined as given 
in eq. ()2.5p . For the Mandelstam variables one finds 

s + t + u = kl + kl+p\+pl=m\ + ml . (4.1) 

We also introduce the following abbreviations: 

h = {k2 - P2f - , ui = {ki - p2f - , (4.2) 

tg = {k2 - P2f - mj , Ug = {ki - p2f - . (4.3) 

When calculating the squared matrix elements of a process with external gluons, 
one needs to sum over all polarization states. In the axial gauge the polarization 
sum for external gluons is given by 

^ (niki) (jiikiY 

with n^^exfi = 0, where is an arbitrary four-vector. This polarization sum as 
well as the polarization vectors obey the transversality relations 

h^^Pr = ni^Pr = (4.5) 

and 

ki^ef^ik^) = eT^{ki)k'; = (4.6) 
for all transverse polarization vectors. This is also true for the polarization sum 



of a massive gauge boson, given in eq. (I3.58p . 



4.2 Squark and Gluino Production Cross Sections in 
SUSY 

For the computation of the partonic cross sections of all processes given in fig. 14.11 
we used the Feynman rules from appendix lA.3[ A way for a simplified treatment of 
external gluons, called ghost subtraction, is treated in appendix lA.2[ The resulting 
squared matrix elements agree with [321] and are given by: 



t-\ ui — mis 



tiui - {mjj - ml)s 



^ufgt NCp ' \ + Agt NCp \/ ^ 

.4^,^ tiui - (m? - m?)s 



^"^'^ Sta 



+ (1 - <5y-) 



/2 



Y.\M''\\gg ^ q-q) = ^njgt 



Co 1-2 



-C. 



K 



1 -2 



srn? 
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J^lM^'l^q^q, ^ qq) = (5, 



+ 4gtmls{NCF 



1 1 

^9 "ff 



2Cf 



tgUg 



+ (!-%) 



iiwi - (m? - m?)s 



i2 



2m^s + tl + ul 




^9t 



SgtNCoil-^-^ 



tg Ug 



2 + 4 



mis 

tgUg 



rn-s 

tgUg 




with iV = 3, Co = N{N^ - 1) = 24, (7^ = {N^ - l)/iV = 8/3 and Cp = 
{N'^ — 1)/{2N). The QCD gauge couphng gs is identical to the Yukawa couphng 
gs- One has to be very careful with the formula given for Yl \-^^\'^{QiQj ~^ QQ)- 
Since the first term derives from the squark production of particles with different 
chiralities the second one comes from the production of equal chiralities. Therefore 
the second term was multiplied by a factor of ^, relative to the first term. As usual 
this symmetry factor has to be taken into account when the partonic cross sections 
are calculated and integration over two identical outgoing particles is performed. 

All matrix elements where calculated completely by hand or with FeynCalc [j^ ]. 
All results agree with the results of [32] and [3^]. In the latter results are also 
given for a hypothetical "polarized LHC". This seems, at least in principle, to 
be an interesting possibility for distinguishing between different models of new 
physics by using spin asymmetries. 

In order to find the leading order cross section at the partonic level we integrate 
the differential cross section 



dtdu 



1 



9{[t-pl][u-pl]-pls)9{s-Am' 



S{s + t + u-pl-p^2)Yl l^^l^ 



(4.7) 



where m = {\/pi + Y^pf)/2. Here Fij is introduced for averaging over the initial- 
states colors and spins: 



QQ 



99 



4(iV2 -1)2^ 



F, 



Q9 



m{N^ - 1) ' 



(4.8) 
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with = 3 for the SUSY-QCD group SU(N). Integration over the two Man- 
delstam variables t and u yields the partonic cross section at leading-order. The 
integration over one of the two variables, e.g. u without loss of generality, is triv- 
ial since integration is performed over the (5-distribution in the double differential 
cross section. 

Integration over t can be performed analytically. Limits arise from the first 
^-function and are given by 



s + mf, 



2 



1 

± - 
2 



(4.9) 



for the case of one outgoing squark and one outgoing gluino in fig. 14.11 f ) . The 
second ^-function is the constraint from the production threshold. It is implicit 
in the following cross section formulae. 

A factor of ^ for identical particles in the final state is included. The resulting 
partonic cross sections agree with [32] and are given by: 



(qiij qq) = <5j^ 
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with 



Li = hi 



s + 2m^ - s(3q 
s + 2w?_ + Sl3q 



In 



s — 2m^ — sf3g 
s — 2m?_ + sfig 



= In- 



In 



s + m 



4m| 



/3a 



4.3 Hadronic Transverse Momentum and Rapidity Dis- 
tributions 

Up to now we have only investigated the colhsions of quarks and gluons. Since 
the LHC is a hadron colhder, one has to take into account the inner structure of 
the proton. The necessary PDFs were already introduced in section [2^21 For our 
analysis of the resulting hadronic cross sections, we choose the PDFs from the 
CTEQ 6 collaboration [9\. 

When hadronic cross sections are calculated, all processes from different initial 
states and with equal final states have to be summed up. 

Hadronic cross sections are calculated numerically since the PDFs can not be 
derived from first principles within perturbation theory and have to be extracted 
from experimental data. Therefore we use Vegas [IB], a Monte Carlo integration 
routine. 

In this chapter we will explicitly calculate the total cross section and the dif- 
ferential cross sections in transverse momentum and rapidity for the process of 
quark-gluon collision producing a squark-gluino pair. Kinematics are very im- 
portant for the comparison of different models of beyond-standard-model physics. 
Since the transverse momentum is invariant under Lorentz boosts along the beam 
axis and can be calculated and measured easily it is a very convenient quantity, 
especially for hadronic processes where the center of momentum is not fixed. 



In analogy to [Sj], we derived a special parametrization for the hadronic inte- 
gration of transverse momentum and rapidity. As an example let us investigate 
the process 

hiiKi) + h2iK2) q{pi) + ~g{p2) . (4.10) 

Ki and K2 are the momenta of the incoming protons while pi and p2 are the 
momenta of the outgoing squark and gluino. Of course we only have to cover 
the case of massless incoming particles. Both outgoing particles are massive. 
Analogously to the partonic Mandelstam variables the hadronic ones are defined 
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by 

Tg = {K2 - P2f -mj, Ti = {K2 - P2? - m? , (4.11) 
Ug = {Ki-p2f -ml, Ui = {Ki-p2? -ml- 

So obviously one can write 

S = X1X2S , tg = X2Tg , Ug = X lU g . (4-12) 

The definitions for pt and y of the gluino are given by 



p't = -^-ml='-l^-ml 2/ = ^ln(ff). (4.13) 



The first relation can be easily motivated by geometrical considerations. The 
formula for y can be shown to be equivalent to 



for massless incoming particles. The second term corresponds to the usual defini- 
tion of rapidity in the partonic frame 



2 VP2 - P2,L 



^P-. = ( ) , (4.15) 



while the first term takes into account that the rapidity in eq. (j4.13|) is given 
in the hadronic frame. Depending on the fractions of momenta xi and X2, one 
has to boost the rapidity for calculating it in the lab frame, since the center of 
momentum is different for every pair of incoming momenta. 

The double differential hadronic cross section in the hadronic Mandelstam vari- 
ables is given by the convolution of the PDFs with the partonic double differential 
cross section 



dTgdUg 



iS,Tg,Ug,Q^) 



f dx, f da:2Xl/^(xl,Q^)x./f (X2,Q^) ^''^'^'^!:'^;"^'^'^ , (4.16) 



where the partonic cross section is given by eq. ()4.7p and multiplication of xi and 
X2 is due to the Jacobi determinant. Distributions in transverse momentum pt or 
rapidity y are gained by using 

dpt dy ^* dTg dUg ^ ^ 

and then integrating out pt or y 

a{S,Q')^ dpt dy-^^ dy dpt . (4.18) 
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The hmits for the integrations for the three variables, xi,X2,y or xi,X2, Pt, in 
order to obtain a differential distribution in pt or y, are simply calculated by using 
the 9 functions and the momentum conserving 5-distribution in eq. ()4.7p . One 
finds 



-xiUg-ml + ml ^__-Tg-ml + ml ^^^^^ 



xiS + T„ ' ' S + U, 



9 



where comes from integrating over the 5 -distribution while the lower limit for 
xi is obtained by setting X2 = 1- The limits for pt and y integration are easily 
derived to be 



pr"(y) = -^J—J(s + ml-miy -Ampcosh^y, (4.20) 



coshy 

y>^^-(p,) = arccosh ( — , ' | . (4.21 




Choosing another parametrization to calculate the cross section is of course 
possible and was numerically done as a check. Since integrations with Vegas are 
performed automatically in an interval from to 1, one has to make a substitution 
for the limits of integration and then multiply the integrand with the appropri- 
ate Jacobi determinant. In cases where integrations are more involved the so 
called histogram method is used, since then one does not have to derive a special 
parametrization. In the later calculation of 2 ^ 5 processes we calculate the dif- 
ferential cross sections for every observable from the four vectors of the particles, 
using any parametrization and within only one run of the program. 

For the hadronic differential cross sections of gluino-squark pair production by 
quark-gluon collision one obtains the distributions shown in figs. 14.21 to 14. 4i Cal- 
culations were performed using the partonic cross sections from section 14.21 and 



using the helicity amplitude generators SMadgraph [Sa], [33] and HELAS [3^ 
These results agree within the numerical accuracy of our integration. The calcu- 
lation of the hadronic cross sections was also performed for the other SUSY-QCD 
processes given in fig. I4.1l a)-e). But since they are not primarily important for us, 
and show no significant new results compared to the squark-gluino production, we 
omit their presentation. For all plots we summed over all quark fiavors except for 
the top-quark. 

As given in the SPSla scenario, we assumed masses to be m^^ = 562.260 GeV, 
= 545.890 GeV and rrig = 606.105 GeV. For the strong coupling we chose a 
value of as(580 GeV) = 0.100375. Masses of all squark flavors are assumed to be 
equal. 

Of course the transverse momentum of squark and gluino are identical for this 
2^2 process. As one can see from [s^], where calculations for the squark and 
gluino production processes are also performed in next-to-leading-order (NLO), 
the NLO-corrections to pt and y distributions are very small. The normalized 
NLO-distribution of transverse momentum and rapidity are quite well described 
by the lowest-order approximation. 
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Figure 4.2: Transverse momentum distribution of gluino and squark in the SPS la 
scenario. 



0.5 



0.4 



0.3 



0.2 



0.1 








-4.0 -3.0 -2.0 -1.0 0.0 1.0 2.0 3.0 4.0 

y 

Figure 4.3: Rapidity distribution of squark and gluino in the SPSla scenario. 
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Figure 4.4: The invariant mass distribution of the squark-gluino pair in the SPSla 
scenario. 
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Therefore is seems legitimate to only use the LO distributions to compare the 
kinematics of SUSY to other models. 

If the PDFs are set equal to one, the shapes of both rapidity distributions 
of squark and gluino are equal. This implies, that the shape of the rapidity is 
governed by the PDFs. In fig. 14.51 we show that significant differences for the 
distributions of particles and antiparticles exist. The red and the green line derive 
from quark-gluon collisions while the blue and the pink line derive from antiquark- 
gluon collisions. The differences are due to the fact that the parton distribution 
function of the quarks is very different from the gluon PDF. This is different for the 
antiquark PDF, which is more similar to the gluon PDF, cf. fig. 12. II This results in 
a large boost for the u- or d-quark and gluon system, because u- and d-quarks are 
more often included in the proton with larger values of x than their antiparticles 
are. This can be seen immediately from the dent in their PDFs. Therefore the 
rapidity distribution for the squarks and gluinos in the case of incoming quarks is 
broadened, compared to the case of incoming antiquarks. The shape of the rapidity 
distributions for squark, antisquark and gluino are differently influenced by the 
PDFs, since in t- and u-channel diagrams the outgoing squark directly couples 
to the incoming quark and the outgoing gluino couples to the incoming gluon. 
As a consequence the rapidity distribution of the antisquark in fig. 14.51 obtains a 
minimum at y = 0. The rapidity distributions of a LHC process always have to be 
symmetric, since the colliding partons can come from both of the protons. There 
is no favored direction at the LHC even for the more complicated 2^5 process. 

The numerical results for the cross sections are given in table 14.11 separately 
for all quark flavors. After summing over all quark flavors, we obtain a total cross 
section of 21.64 pb. 
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Figure 4.5: Rapidity of squark and gluino from quark-gluon collision and anti- 
squark and gluino from antiquark-gluon collision in the SPSla scenario. 



Chapter 4 



SUSY- and UED-QCD Processes 



39 



4.4 The Production of a Kaluza-Klein Gluon-Quark 
Pair in UED 



In this section we compute the 2^2 cross section for KK-gluon and KK-quark 
i9*,Qn,i) production. The diagrams taken into account are presented in fig. 14.61 
First calculations in the literature were performed in 23] and Their results 
were corrected in [4]- Both calculations were done for equal masses of KK-gluon 
and KK-quark. This is reasonable for a UED calculation since masses are nearly 
degenerate in a UED spectrum. But since we want to compare SUSY and UED for 
the SPS la mass spectrum, we derive a more general cross section with different 
masses of the final state particles. As a check we could show that in the limit 
of equal masses for the KK-quark and KK-gluon, the result of 0] agrees with 
our calculation. Using the Feynman rules from the appendix and the polarization 
sum —gfj,i, for gluons and eq. (|3.58|) for KK-gluons, we obtain the squared matrix 
element given by 

i:\M{qg^qn,ig:)\'- ^'^^^ 



1152 s tl, u 



X (4tiq.ml + {mL{4tg 



+ng* + s^iiq.jm^, + TTig* tg*{—2m^* — 2uq*mq, + tg*Uq*} 



2uq*} + 4tg*Uq*}mg* + 2{smq* — tg*Uq*mq* 

(4.22) 

^, and tg = t — m?g* . The masses m^* and nig* here denote the 
masses of the KK-quark and KK-gluon. The squared matrix element in the limit 
of equal masses of the outgoing KK-gluon and KK-quark, as it is the case in a 
typical UED mass spectrum, is given by 



with Uq = u 



J^\M{qg^ qn,ig*n) 



-1 



5s" 
12t'2 



llsu' 5ti'2 w 



/3 



+ 



with u' = u — and t' = t 
eq. (I4.22p . according to eq. ()4.7p . and obtain 

a^(Q)vr 



o- {qg qn,ign) 



36 m^, 



16/umg. -|- 4{3 (/i — 3u)m'^* + r]}m't* 



+ 2 {-(9 {fi-iy)ml, +ri)ml, +115^ + 25^ {2fi + 9i^)}ml, 



+ 365^^ - 2ml ?7 + m^,s {3^ + 4s/x} - 2ml {5/i - 9i^} 



(4.23) 



with 



m^ + {rriq — s}'^ — 2m^{mg* + s} 



In [nigt 



ni„ 



In ( m^* 



ni„ 



In (m^* 



niq* + s 







In ( m^* 



r] = 16C- is/j, + 9si^ . (4.24) 
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Figure 4.6: LO contributions to the production of heavy KK-quarks and heavy 
KK-gluons in UED. 

Calculating the kinematic distributions in UED for the parameters of the SUSY 
SPS la benchmark point using Madgraph, one obtains the graphs given in figs. 14.71 
to 14.91 We checked the agreement with the given formulae numerically. The 
transverse momentum of KK-quark and KK-gluon are equal due to momentum 
conservation, as it is also the case in our SUSY calculation. The invariant masses 
of the outgoing particles in SUSY and UED show a different behavior. This is 
shortly discussed in section IT5l The difference between the squark and KK-quark 
rapidity distributions are discussed in section 14.61 in more detail. 

Later a comparison of the hadronic SUSY and UED distributions will be per- 
formed in chapter [6l The numerical results for the cross sections are given sepa- 
rately for all quark flavors in table 14. 1[ From these values one can see that UED 
cross sections are much larger than SUSY cross sections. Accordingly, the total 
cross section itself could be used to discriminate UED and SUSY at the LHC, 
provided that the masses of all participating particles are known. After summing 
over all quark flavors, we obtain a total cross section of 170.08 pb. 
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Figure 4.7: Transverse momentum distribution of KK-gluon and KK-quark with 
masses of the SPSla scenario. 
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Figure 4.8: Rapidity distribution of the KK-quark and KK-gluon with masses of 
the SPSla scenario. 
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Figure 4.9: The invariant mass distribution of the KK-quark-KK-gluon pair. 
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Table 4.1: Cross sections for all quark flavors for quark-gluon collision in SUSY 
and UED 
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Figure 4.10: Comparison of the hadronic invariant mass distributions of the given 
particles in SUSY and UED. 



4.5 A SUSY-UED Comparison at the Threshold 

Investigating the behavior of the cross section close to the threshold is interesting 
due to its strong energy dependence. It can be expected that SUSY and UED 
models differ significantly in this region. As we can see in fig. 14.101 where the 
invariant masses of a squark-gluino and a KK-quark-KK-gluon pair are shown, 
differences right behind the production threshold are obvious. This could be due 
to a different threshold behavior of the partonic cross section in both scenarios. 
Therefore it should already be visible at the partonic level. 
As in [i^], the so called phase-space suppression factor /3, is defined as 



/ ps^^, (4.25) 

y s-[mq-mgY 

Close to the threshold, i.e. when the produced particles have a small velocity in 
their center of mass system, /3 ^ 1, one can give analytical expressions for the 
expansion of the cross section in (3. The expansion of the SUSY cross section for 
quark-gluon collisions, as given in section W?I[ yields 



^susT°''(^5 Q9) = a'iiQ)- — 



(4.26) 



This agrees with the result given in [33|- In the UED scenario one obtains for the 
equivalent process, expanding eq. ()4.23p 



B approx 



{qg ^ qn,ig*n) = (4.27) 

, {^ni^ + 34mgm| + JAm'^m'^ + 81m?mg + 45m|) 

= a,(Q)7r/3 ^ ^ ^ J — ^ ^ . 

9mg[mg + mg)^ 
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Figure 4.11: Partonic leading-order cross section in UED for production of a KK- 
quark-KK-gluon pair. 
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Figure 4.12: Partonic leading-order cross section in SUSY for production of a 
squark-gluino pair. 
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Figure 4.13: Comparison of the UED and SUSY cross section threshold expan- 
sions. 
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In fig. l4.11l and fig. l4.12l we compare the partonic SUSY and UED threshold ex- 
pansions to the full partonic cross sections. In fig. 14.131 we compare both threshold 
expansions. The masses of the outgoing particles are given by rUg = 606.105 GeV 
and niq = 562.260 GeV. The factorization scale is set to Q = 580 GeV. The 
quantity r] is defined by 

V = 7 V? - 1 • (4-28) 

[mq + rugy 

There are obvious differences between the partonic SUSY and UED cross sec- 
tions. While the SUSY cross section peaks behind the threshold and decreases 
for s — > oo, the UED cross section does not decrease. This is in contradiction to 
unitarity. But since we did not take into account higher KK-excitations, due to 
suppression by their higher mass, the result is only meaningful below the second 
KK-excitation, i.e. not too far above the threshold. Therefore it seems to be rea- 
sonable to compare UED and SUSY at the threshold. Directly finding the second 
excitation level at the LHC would, of course, be a clear sign for the theory of 
Universal Extra Dimensions being realized in nature. As one can see easily, the 
UED cross section increases much faster than the SUSY cross section. Of course 
this partonic quantity can not be measured at the LHC. But it shows that typical 
UED cross sections are usually much higher than typical SUSY cross sections, 
assuming the same mass spectrum. 



4.6 A SUSY-UED Comparison of Angular Distribu- 
tions 

As we find in fig. 14.31 and fig. 14.81 the rapidity distributions in SUSY and UED 
are quite different. These hadronic rapidity distributions are all shown again in 
fig. l4.14l One finds that the gluino and KK-gluon rapidities are quite similar, while 
there are differences for the outgoing squark and KK-quark. For the hadronic cross 
section the sum is taken over all quark fiavors u, d, c, s and b. 
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Figure 4.14: Comparison of the hadronic rapidity distributions of the given parti- 
cles in SUSY and UED. 
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A plausible explanation is found by analyzing the angular distributions of the 
partonic cross sections in the CMS, given in fig. 14.151 and fig. 14.161 For the case of 
a moderate center of mass energy, e.g. ^/s = 2.0 TeV, one finds that the angular 
distributions of the squark and KK-quark in SUSY and UED have a very similar 
shape. For a high center of mass energy, e.g. ^/s = 8.0 TeV, one finds increas- 
ing differences between SUSY and UED. These differences could be due to the 
momentum dependent three gluon vertex in UED, while the corresponding gluon- 
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Figure 4.15: Distribution of the angle between the squark/ KK-quark and the 
beam axis in the CMS for ^/s = 2.0 TeV. 
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Figure 4.16: Distribution of the angle between the squark/ KK-quark and the 
beam axis in the CMS for ^/s = 8.0 TeV. 



Chapter 4 



SUSY- and UED-QCD Processes 



47 



gluino-gluino-vertex in SUSY is independent of the incoming momenta. Instead 
the squark-squark-gluon vertex in SUSY is dependent on the squark momenta. 
Therefore it would be very surprising if the angular distributions in SUSY and 
UED had an equal shape for all center-of-momentum energies. 

The distributions are computed using the relations for the partonic cross sec- 
tions in chapter [2j The angle 9cm,i between the beam axis and the outgoing 
squark, respectively the KK-quark, is given in the center-of-momentum frame as 
defined in eq. ()2.7|) . In the CMS, the outgoing gluino, respectively the KK-gluon, 
has an angle of 9cm, 2 = 180° — 6cM,i- Since the hadronic rapidity distributions 
include this information from the angular distributions on parton level, it is plau- 
sible that the rapidity distributions can differ on the hadronic level. The u- and 
(i-quark are comparatively often included in the proton with a high momentum 
fraction. This yields a high center of momentum energy s, where angular distri- 
butions in SUSY and UED are somewhat different, which can lead to differences 
in the hadronic rapidity distributions. 

As we will see in chapter El the difference in the hadronic differential cross 
sections of the rapidity for SUSY and UED yields a difference for the kinematics 
of the decay products of the squark, compared to the decay products of the KK- 
quark. 
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At the LHC one will not observe supersymmetric particles directly. Instead people 
are searching for their decay products. It is important to choose a signature that 
is as clear as possible, considering the background from Standard Model processes. 
Therefore decay chains play an important role in the detection of supersymmetric 
particles. To be certain that one observes SUSY at the LHC one has to measure 
the spins of the decaying particles. In a comparison of the kinematics of an extra- 
dimensional and a supersymmetric theory, differences in kinematic distributions 
should show up, due to different spin correlations in the applied models. In a 
theory with extra dimensions all "extra-dimensional" partners, the Kaluza-Klein 
particle towers, have opposite spin statistics compared to the supersymmetric part- 
ner particles. In our calculations we assume that all masses are already measured. 
Masses can be found by using kinematical distributions like two- and three-particle 



invariant mass distributions [39|, 40(]. In this chapter we discuss the kinematics 



of decay chains and explain the setup of our program. 

5.1 Decay Kinematics and Phase Space Generation 

If a particle with momentum p decays into n particles with momenta pi, . . . ,pn, 
the phase-space element, described by (3n — 4) independent variables, is given by 



ncfpi 

.1=1 



S^'Hp-^Pi). (5.1) 



i=l 

Hence, the phase space for one particle decaying into two particles reads 
d<^{p ,m^,m2) = J ^^^^^ ^P~P^~P2) 

1 2n 1 

X^{p'^,ml,ml 



8p2 



y # y dcosO , (5.2) 



-1 



where A is given by eq. (j2.10p . The momenta of the outgoing particles in fig. 15.11 
are given in the rest frame of the decaying particle by 

\Pi\ = \P2\ = yPi = ■ (5.3) 

* 1\J p^ 
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P2 



Figure 5.1: Our decay chains are built of 1 — > 2 particle decays. 

In order to find the momentum of the particles in the lab frame one has to boost 
them to the frame where the decaying particle has the momentum p. The usual 
Lorentz boost along the z axis, from the rest frame into a frame moving with /J^, 
is given by 

L{7,P,)pi, (5.4) 



with L{'y,Pz) given by 



P'l 



Hi, (3, 



/ 7 

1 

1 

V-7/3. 






7 / 



(5.5) 



and 



= ^ 
m 



m 



(5.6) 



For a chain of multiple decays we need a more general boost matrix into a frame 
moving with (3 in the direction |4| , given in 41] by 



i(7,/?) 



V- 



7 




-7/?x 




-7/3, 




-7/3. 


7/3x 


1 + 




(7- 




(7- 




7/3, 


(7- 




1 + 


(7-l)/3^ 


(7- 


-l)/3j//3z 


7/3. 


(7- 


-l)/3:rA 


(7- 




1 + 


(7-1)/?.' 



(5.7) 



Our phase space is generated by calculating four vectors from transverse momen- 
tum and rapidity. Four vectors simply transform by eq. ()5.7p under arbitrary 
boosts. Transverse momentum and rapidity also have a simple transformation 
behavior under Lorentz boosts, while a parametrization in angles and modulus 
of the momenta of the particles has a very complicated behavior under general 
Lorentz boosts. 

Since a real cascade decay involves more than one decay, we need to extend 
the phase space element towards a particle decaying into n particles, by joining 
numerous two-particle decays. A general discussion on that can be found in a 
chapter on multiparticle production in [rj. 

For the phase space element of a chain as shown in figure [521 one obtains the 
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Figure 5.2: The process Pa + Pb Pi + ■ ■ ■ + Pn as a sequence of two-particle 
decays. 



recursive relation 



/ 



1=1 K. 1=1 



d^Pn 



2En 



Rn-lip-Pn). (5.8) 



When one inserts 

1 = J dMti6{Mti-kl_,) (5.9) 

1 = j d^'kn-lS^^Hp-Pn-kn-l) (5.10) 

into eq. (j5.8p and uses the phase space for the two particle decay eq. ()5.2p . one 
obtains 

(M„-m„)2 

Jd^i^n = I dM^_,R2iklkl_„pl)Rn-i{M^_i) (5.11) 



Mn-1 



(M„-m„) 



f , f A2 (M^,M^ i,m^) , „ , 

J dMti J d^n-i ^ "^ i?^_i(M^,) , 



with the trivial limits 

lXi = mi + ... + mj , M^_i = (p - Pn)^ = (Pi +P2 + ■■■ + Pn-i)^ = kl_i . 

The decay is split into a two particle decay with the momenta pn and /cn-i and the 
phase space for the decays of kn-i into the particles with momenta pi, ■ ■ ■ ,Pn-i- 
Integration over all invariant masses in the intermediate states has to be per- 
formed, since the particles do not necessarily need to be on the mass shell. The 
integration limits are trivial since the decay only takes place if M„ > Afn-i + 
and M„_i > Hn-i- For the case of n = 3, the corresponding phase space element 
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exphcitly reads 



{M2-mzY ^ 



X2{p'^,pl,kl) xHklpl,pl) 



8M| 8M| ^ (5.12) 

27r 1 27r 1 



^1 



(icos0i y d</)2 J dcos@2- 

0-1 0-1 

In our apphcation the decay chain follows the production of the decaying particle 
in a 2 —s- 2 scattering process. Therefore the decaying particle's momentum itself 
is not on-shell and also has to be integrated over. 



5.2 Matrix Element 

The matrix element for a decay chain, in principle, is nothing special. But of 
course it is hard to square, since formulae are quite long even for a leading or- 
der calculation. Therefore we explain in this chapter how the matrix element is 
calculated numerically. We also describe what is done to avoid divergences from 
propagating particles. 

In order to keep the perturbative error to a tolerable level, one takes some 
effects of higher orders of perturbation theory into account. One example is the 
running of the coupling which was already mentioned. Since usual propagators 
of massive particles are divergent at the pole, one has to take into account the 
width of the particles propagating in the graph. Otherwise divergences in the 
s-channel diagrams occur. Due to the perturbative expansion 



+ ... 



one finds for the propagator 



+ 



p2 — niQ p^ — rrin 



p^ — rriQ — M^{p'^) 



' 9 2 
— ITIq 



(5.13) 



When the propagating particle is unstable, M'^{p'^) acquires an imaginary part 
and can be written as 



iZ 



p^ -ml- Re M'^{p^)-iZ Im [p^] 



(5.14) 



where Z is called the field renormalization constant. Therefore, close to the pole, 
i.e. when the particle is nearly on-shell, one finds the following dependence of the 
cross section 

1 



a oc 



p2 _ _|_ ^ jyj p 



(5.15) 
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with = ttiq + Re M'^{m^) and 

Z 



r = ImM\m^) . (5.16) 



The propagator including the physical mass m and the particle width F is called 
Breit-Wigner propagator. A more detailed discussion on how to treat divergences 
by regularization can be found in [3]. 

If we plot the value of eq. (|5.15|) against the squared momentum of the particle, 
where the ratio of the particles mass m divided by its width T is small, one can 
see that the propagator strongly peaks at the pole. Therefore factorizing the ma- 
trix element at the Breit-Wigner propagator and forcing the propagating particle 
on-shell seems practicable to make computations easier. This approximation is 
also called "Narrow Width Approximation" (NWA). It corresponds to the multi- 
plication of the cross section for on-shell production of the decaying particle, with 
the branching ratio of the decay. Integrating over eq. (I5.15p . using the residue 
theorem, yields 

— -|- zl ml 
if r <C m. Since T for the decay of a particle p into the final state f is defined as 

1 1 



I d^,^n\M{p ^ ft, (5.18) 



(27r)3'^-4 2mi ^ 
one finds for the example of a 1 — > 3 process 

X {2lr 2m\, Ef I d<i>k2^p,,p;\M{k2 ^ P2,P3)\^ (5.19) 

= T{ks ^ pi,k2) ■ Br{k2 ^ P2,P3) , 



with momenta defined as in fig. (j5.2p . Applying this NWA is problematic if the 
mass of the decay products is close to the mass of the decaying particle. This is 



due to threshold effects and was investigated in 42] and [43 1. 

In our calculations we do not assume that particles are totally on-shell. But 
making use of the fact that the matrix element far from the pole is small, we 
only have to integrate over a certain width around the pole. We will later give a 
numerical justification for this approximation. 

This Breit-Wigner propagator only has to be taken into account in s-channels, 
where particles can be on-shell. Therefore the s-channel is effectively lowered in 
its order in the coupling constant Os by division by F from eq. ()5.17p . since it 
is proportional to the squared coupling constant. In decay cascades involving 
t- and u-channels, the squared momentum is far from the pole. Therefore the 
contribution of such diagrams is not dominated by the width of the Breit-Wigner 
propagator. In this sense these contributions to our decay can be regarded as a 
NLO contribution and can therefore be neglected in our LO calculation. This is 
discussed in detail in section I 



Since we do not want to square the matrix element of a cascade decay process 



by hand, we use the helicity amplitude formalism 4J]. For our calculations we use 
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the tool SMadgraph [3a] , [37[ , which generates a Fortran code calling HELAS [3^ 
a helicity amplitude generator. Madgraph evaluates all possible helicity combina- 
tions for all topologies for the given external particles and produces a Fortran file 
with HELAS calls. HELAS calculates and squares the matrix element numerically, 
using Dirac four-spinors given by 



u\p) 



^x{p) Xxif) 



with 

UJ: 



Vet\p\ . 

Here the two component helicity eigenstates are given by 



x+ip) 



\P\ +Pz 

^/mM+pz) yp- + ^py 



X-{p) = ( + '^M (5.21) 



for \p\ / —pz- For \p\ = —pz one uses 

x+{p) 



X-{P) = ( • (5.22) 



These helicity eigenstates satisfy 

(T • p 
IK 



Xxip) = \Xx{p) (5.23) 



with A = ±1. 

Checking the calculations of Madgraph and HELAS is easy since one can have 
a direct look at all vertices. HELAS performs calculations in a model independent 
way. Since the structure of a vertex is given by the spin and polarization properties 
of the particles coupling to it, HELAS provides general vertices for scalars, vectors 
and fermions. We want to give a short example of a helicity amplitude of a 
2 — > 4 matrix element. In fig. 15.31 one of the diagrams for the gluon-gluon to 
sbottom-anti-bottom, sbottom-anti-bottom process is presented, using the HELAS 
abbreviations for the vertices, which can be found in [38.] . The helicity matrix 
element is given by 



2 2 J -m~+ imgY-g 

,. 1-75 , . l + 75\ i(^2 + "^) (■ 1-75 , . 1+75 
^ 2 2 J q^ — + tnigYg \ 2 2 

i{i^+m) 1-75 , . l+75\ A, 



- m| + im-gYg 



m) m^r^ + i92—:r^ <(m) • (5.24) 
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Figure 5.3: One of the diagrams contributing to the process gg bbibbi 

Here A and A' denote hehcity eigenstates, according to eq. (I5.20p . The matrix 
element is evaluated numerically for every possible combination of helicities and 
polarizations of the external particles. 

In order to calculate the cross sections we need all couplings and masses. In 
the formula above, the coupling gi is the coupling to the chirality left fermions 
while the coupling g2 belongs to the chirality right fermions. They are simply 
derived from the strong coupling constants as given by the Feynman rules. 

As we mentioned earlier, we assume the masses of the SUSY particles to be 
already known. In our calculations we use SDECAY [i^ which calculates branch- 
ing ratios for the supersymmetric particles, making use of the spectrum generator 



SuSpect 46|. SDECAY puts masses, mixing angles and branching ratios in a 



SLHA format file 47[, which is read into our program and used by our matrix 
element and phase space generator. In SuSpect a phenomenological MSSM with 
22 free parameters is implemented. When parameters are set by the user, renor- 
malization group running is performed, going iteratively up to the GUT scale and 
down to the electroweak scale, taking into account threshold effects of all particles. 
All masses and couplings are then derived for the electroweak scale, also including 
higher order effects. 

Our whole program is build up step by step on the basis of the first version 
of the simplest 2 — > 2 scattering cross section program. But since the number of 
degrees of freedom (3n — 4) increases with the number of the outgoing particles 
n, the integration gets more complicated. Therefore it is mandatory to introduce 
mappings, forcing the integration routine to evaluate more phase space points in 
the region where the value of the integrand is large. The integration variable x is 
mapped to a set of random numbers y by 

x = h{y), Os^y^l. (5.25) 

The integral can then be calculated as 

/ = / fix)dx = C nhiy))^dy = t ^dy , (5.26) 



dy Jq g{h{y)) 

where g{h{y)) is called the density. The integral is then calculated by sampling 
the integral N times and averaging as 
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By mapping the integration variables x in a way that is smoother than 

f{x), convergence of the integral is improved and more points are evaluated in 
those regions where f{x) is steep. 

The most important mapping in our program is the one for the Breit-Wigner 
propagator , as given in by 

h{y, - imT, x^in, Xraax) = rriT tan [zi + (z2 - zi)y] + . (5.28) 
This results in a density given by 

mV 

g{y, va^ - imF, Xmin, a^max) = 7 777 n7n-, (^-29) 

[Z2 — zi) [[X — m'^y + m^i 

with 

( ^ rain /max ~ \ /c o^^ 

Zij2 = arctan I — 1 . (o.oUj 

Leaving out this mapping results in a very inefficient and slowly converging inte- 
gration. The function we integrate over, the squared Breit-Wigner propagator, is 
given by 

f{p) = -7—2 fF2 ■ (^•^^) 

As we will see in section 15. 3[ our limits for the momentum of the propagating 
particle are determined by its width 

— n-r + m ^ -p ^m + n- F, (5.32) 

which yields 

f (m± nF)^ — 

Zi/2 = arctan ( I . (5.33) 

Calculating the integrand in eq. ()5.26p one effectively finds a constant function 

fiHy)) Z2 - zi 



g{h{y)) mV 



(5.34) 



which has to be integrated from zero to one. This constant function can be 

integrated easily while the squared Breit-Wigner propagator is strongly peaked at 

2 2 
p = m . 

Other mappings, e.g. a mapping for the integration over pt using pt = ^, are 
also included. Angles are mapped as 

cos6' = 2y — 1 and cj) = 2-Ky . (5.35) 



Since our Monte Carlo routine always integrates from to 1, Jacobians have to 
be included for different integration limits. 
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5.3 Testing of a Decay Chain Program 

In the following we explain the testing of our phase space generator. Beginning 
with a usual parametrization of a 2 — > 2 process in terms of y and pt and inte- 
grations over xi and X2, we extend the phase space by attaching a 1 ^ 2 decay 
to one of the final state particles. Then the squared momentum of this decaying 
particle has to be integrated over. We make use of the fact that all intermedi- 
ate state particles have a small width compared to their mass, when we discuss 
the emerging topologies. For a three particle final state there are three different 
topologies: 




b) 



Pi 

P2 

P3 



In our calculations we take into account the topologies of b) and c). Since the 
momentum in the t-channel is far from the pole, the topology a) is not enhanced 
by the Breit-Wigner and therefore suppressed, compared to the other topologies. 
The following topologies can be found for four particle final states: 



a) 



c) 



k2_ 



Pi 

P2 
P4 



Pi 
Pi 
P3 
Pi 

Pi 
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The topologies a), b) and c) are again neglected since their momenta are far from 
the poles. In general, all other topologies are taken into account. Depending on 
the mass hierachy of the intermediate particles, it is possible that topology g) is 
also suppressed. This happens when an intermediate quasi-on-shell particle, i.e. a 
particle with a momentum within the limits of eq. (j5.32p , is followed in the decay 
cascade by a heavier quasi-on-shell particle. It is not possible for both particles 
to be close to their mass shell at the same time. 

The topologies of our special 2 — > 5 process will be investigated in chapter [6l 
Leaving out all effective NLO topologies, we are left with s-channel decays and can 
compare, e.g. the 2^5 process to the 2 — > 4 process multiplied by the branching 
ratio of the last decaying particle, because the decaying particle is nearly on-shell. 
To give a feeling for the contribution of the on-shell momenta, we give some total 
cross sections in table 15.11 Calculations were performed for the 2^2 processes 

u g g and g g ^ g g . (5.36) 

To obtain a 2 — > 5 process, we attach step by step the decays 

(1) g^b, hi (2) b^b,N2 (3) u^u,Ni (5.37) 

to the first process in eq. (I5.36p . This was done for an intermediate ul and u/j. 
We also repeated the procedure for the second process, given in eq. (I5.36p . Since 
the masses in this process are different, it gives us a second check for the generated 
phase space. 

The columns in table 15.11 represent three different integration limits for the 
intermediate momenta. Due to the Breit-Wigner propagator, contributions are 
only large at the pole and therefore we integrate an integer multiple of widths, 
here denoted n ■ T, around the pole, i.e. our explicit integration limits are given by 

{-n-r + mf ^ / ^{m + n-rf. (5.38) 

The diagrams contributing to the cross sections are shown in fig. 15.41 Though other 
diagrams and more complicated topologies in principle exist, their couplings were 
commented out in Madgraph in order to make a comparison with the branching 
ratio possible. Row one to three in table 15.11 correspond to intermediate bi and 
Ul while row four to six correspond to intermediate 6i and ur. If there are more 
constellations for the final state, e.g. the commutation of the b and b jets, we 
left them out since we just want to understand the behavior of our phase space 
here. We only included the bi squark and left out 62 in the intermediate state. 
Therefore these calculations are by no means measurable results, though we could 
easily include all particles. 

The masses for this arbitrarily chosen scenario are given by rrig = 746.01 GeV, 
ruu^ = 692.72 GeV, m^^ = 625.12 GeV, uin, = 93.83 GeV, = 191.15 GeV. 
For the strong coupling we used as = 0.1003750. The used particle widths are 
given by Tg = 8.9834 GeV, Ta^ = 7.7452 GeV and F^^ = 6.3995 GeV. The 
branching ratios corresponding to this spectrum are calculated with SDECAY 
and are given below. 
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Figure 5.4: Diagrams contributing in LO to the 2^5 process. 



Process 


o-„=5 in pb 


o-„=i5 in pb 


o-n=25 in pb 


o-n=30 in pb 


u g ^ ULbbi 


174.0(1) 


182.0(1) 


183.9(1) 


184.5(1) 


u g ^ uib b Ni 


53.13(3) 


55.57(3) 


56.17(3) 


56.34(3) 


u g ^ u Nib b N2 


0.4152(2) 


0.4532(2) 


0.4616(2) 


0.4638(2) 


ug ^URbbi 


186.0(1) 


194.5(1) 


196.6(1) 


197.1(1) 


u g ^ URbb Ni 


56.79(3) 


59.38(3) 


60.02(3) 


60.20(3) 


u g ^ u Nib b N2 


53.13(2) 


58.06(3) 


59.18(3) 


59.49(3) 



Table 5.1: Cross sections for an exemplary calculation. 
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We again stress that our result does not take into account off-shell contributions 
in general, since beyond the 2^2 process we only consider s-channel decay 
topologies, while we still allow some deviation of the squared momentum from the 
pole. In a region not too far from the pole, as given by eq. ()5.38p . contributions 
to the total cross section can be quite significant. From table EH] we find, that the 
main contributions to the process come from intermediate on-shell momenta. In 
our later calculations we integrate out a region with n = 25, since then only small 
further changes in the cross sections appear. We also find that our integration is 
equally stable for all cases. Integration errors are given in brackets for the last 
digit. 

For checking the phase space, we use the narrow-width-approximation. For 
the case of n = 5 and ul, taking the branching ratios from SDECAY, we multiply 
the 2^3 result by 

BR{bi N2b) = 0.305532 (5.39) 

to obtain 53.17 pb, which is close to the calculated value of 53.13 pb. For the 
decay of the ul squark we multiply by 

BR{uL ^ N2b) = 0.0082336 (5.40) 

and obtain 0.4375 pb, compared to the calculated cross section of 0.4152 pb. 

For the case of n = 25 and ur, we multiply the 2 — > 3 result by BR{bi — > b) 
to obtain 60.06 pb, which is close to the calculated value of 60.02 pb. For the decay 
of the ur squark we multiply by 

BR{uR N2b) = 0.995188 (5.41) 

and obtain 59.73 pb, compared to the calculated cross section of 59.18 pb. This 
approximation is quite good for larger and smaller n. Obviously ur couples much 
stronger to the neutralino N2, which is due to the mixing angles of neutralinos 
and squarks in the given scenario. 

As a last useful check for the extensions of our program, we used various plots 
to qualitatively understand the behavior of the program and check for consistency. 
The shape of transverse momentum and rapidity distributions as well as invariant 
mass distributions have to stay nearly unchanged for the intermediate particles 
when the phase space is enlarged by an additional 1^2 process, since only small 
off-shell infiuences are included. This was tested and various distributions will be 
presented in the next chapter. All these agreements confirm that our phase space 
generator works correctly. Since for the case of UED calculations we only have to 
exchange the matrix element, phase space checks were only performed for SUSY 
calculations. 
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Chapter 6 

A SUSY-UED Decay Chain 
Comparison 

Finding like-sign dilepton signatures is a promising strategy to confirm new physics, 
as explained in [s^. Like-sign signatures are expected to emerge from processes 
involving Majorana fermions, e.g. the gluino in SUSY. Processes like qq/gg — > gg 
and their decays are candidates for this discovery. The Majorana gluinos decay to 
qq or qq, leaving like-sign leptons in the final state. But the idea of using like-sign 
dileptons to claim the existence of a heavy Majorana fermion is incomplete. If the 
particle decaying in the gluino-like cascade is a boson in the adjoint representation, 
e.g. a KK-gluon in UED, like-sign signatures also arise. 

To solve this problem and show that SUSY is the theory realized by nature, 
one has to show that the decaying particle is really a fermion. This can be done by 
comparing two scenarios of different spin assignments for the gluino-like particle. 
UED is a popular candidate to compare to SUSY. It has exactly the same spin 
assignments as the Standard Model for all partner particles. Since the UED mass 
spectrum is typically very different from a SUSY spectrum, there are many ways 
to distinguish between UED and SUSY. The mass spectrum itself, the discovery of 
higher KK-excitations, ratios of branching fractions, threshold behavior and cross 
sections can be expected to be very different for both theories. Since we want to 
use UED as a toy model with different spin assignments for verifying SUSY, we 
are not interested in the total cross section of a UED process. As usual we extract 
spin information from angular correlations. Therefore we divide all distributions 
by the total cross section and exclude influences from different masses and coupling 
constants. An entanglement of the spin information with the couplings of the left 
and right handed sfermions in our decay chain can not be excluded. 

We show in this chapter that the spin information can be extracted from 
kinematic distributions. Differences in the boost distributions and angular cor- 
relations are studied. We also investigate the origin of the differences in angular 
correlations. 

Various recent publications deal with the issue of comparing new-physics sce- 
narios to find the spin of new particles emerging at the LHC. Knowing spin and 
masses is crucial for knowing the Lagrangian. While masses can be extracted 
from various edges and thresholds of invariant mass distributions [s^, [iO], the 
measurement of the spin is non-trivial. 
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The idea first came up in [48|] , where it is demonstrated that spin information 
can be extracted from angular distributions and invariant masses. This idea is 
used in [2] to perform a comparison of UED and SUSY decay chains to determine 
the spin of the squark. Analytical results for the invariant mass in the decay of a 
squark are given. The leading order UED-QCD cross sections are also presented 
in this paper. 

Due to the like-sign dilepton argument the spin measurement of the gluino is 
especially interesting. Therefore LHC decay chains involving a gluino and KK- 
gluon are analyzed in [H] , in addition to our decay chain also including the further 
decay of the NLSP into the LSP and two leptons. For the processes of gluon-gluon 
and gluon-quark collisions various asymmetries of outgoing leptons and bottom 
jets are studied. It is shown that by using these asymmetries the spin of the gluino 
can be determined. Significant differences in the angle between two outgoing 
bottom jets are found. It is claimed that this difference is mainly due to the 
different boost of the gluino and the KK-gluon. 

In a recent publication the question arises if the boost of the gluino-like 
particle or the different helicity structure of the couplings is responsible for the 
kinematics in the final state of a decay chain. This is studied for the three-particle 
decay g ^ b b Ni. It is argued that none of both influences can be excluded from 
the calculations. Even in the case that masses and spins are treated as independent 
parameters which are used to fit one model to the other, it is shown that SUSY 
and UED can always be differentiated by invariant mass distributions of the decay 
products. It is explained that the longitudinal modes of the neutralino, if highly 
boosted, are mainly responsible for the characteristic differences. Calculations are 
performed assuming rrig > rriq for all squarks. 

In this chapter we calculate the gluon-quark collision and analyze the differ- 
ences in the gluino and KK-gluon boost and the angle between the outgoing jets. 
We leave out the decay of the NLSP. In the SPS la scenario, the first neutralino 
A^i is the LSP and the second neutralino N2 is the NLSP. Together with these 
neutralinos, two 5-jets and a light quark jet are outgoing. In comparison to [H] 
and [sl, we additionally try to find out numerically if the impact of the boost on 
the angular correlations between the 6-jets is significant. This is done by mapping 
the boost distributions of the SUSY gluino and the UED KK-gluon onto each 
other and thereby taking out the effect of the different boosts. 



6.1 Decay Chains in SUSY and UED 

The decay chain we want to consider is presented in fig. 16.11 Since jets of quarks 
and antiquarks are not distinguishable one has to sum over all flavors of quarks 
and antiquarks. Using SDECAY we obtain different masses for up- and down-like 
flavors and left- and right-chirality of squarks as given in table 16.11 We therefore 
include all possible combinations of ul, ur, and respectively dL, cIr, and 61, 62 
in our calculations. Since the down-like coupling to the LSP is different from the 
up-like one, it is essential to evaluate the matrix elements separately for these two 
kinds of flavors due to its dependence on the weak isospin. This dependence is 
not present in the case of the 2^2 process in chapter 4. There one can simply 
sum over all PDFs and multiply them by the matrix element of the 2 — > 2 process 
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Figure 6.1: Diagrams contributing in LO to the 2 — > 5 process in SUSY. Ah hght 
incoming quark flavors can not be distinguished and have to be summed over. 
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o u o I -pdrxici6 


U H/U-pcLi tlClG 


iTLcLSS in GgV 


WlQIIl in \jrG V 






606.105053 


4.51087328 


11, T 


7/1 1 


562.260527 


5.45839994 






545.890747 


1.14435575 


dj 

'-''Li 


d^ 1 


567.765846 


5.25959596 






545.619560 


0.28594360 


hi 


bi 2 


516.907879 


3.76415860 


In 




546.235967 


0.85335990 


Ni 


51 


97.1850209 


0.00000000 


N2 


A3S 


180.301343 


0.02205761 



Table 6.1: Masses and widths for the SUSY particles used in our calculations. 

if masses of all flavors are assumed to be equal. 

Exchanging the b- and 6-jets in the final state of the 2^5 process is also 
necessary. The sbottoms are then substituted by anti-sbottoms. We do not con- 
sider incoming bottom-quarks since there would be an additional bottom-quark 
in the final state, leading to additional permutations with the other bottoms in 
the final state. This leads to a highly increased number of diagrams that have to 
be evaluated. Moreover this is not necessary from an experimental point of view 
since it is practically possible to differentiate between a bottom jet and a light jet, 
i.e. u,d,c and s-jet. 

Since the same final state can also be reached by exchanging A^i and A'^2i one 
has to exchange them as well. Both particles can in principle also couple to the 
other vertex. But since we do not claim to perform a calculation with a directly 
measurable result, we leave out the exchange of A^i and N2 and respectively the 
KK-particles and ^3,1. This does not harm gauge invariance explicitly since 
there are no gauge parameters involved. Therefore the calculation is a reasonable 
but not directly observable one. As already stated, we are interested in the origin 
of the angular correlations of the bottom jets. This is a theoretical question 
and exchanging these particles would result in doubled run time, not providing 
a much deeper insight into the angular correlations. In the SPS la scenario the 
branching fraction of the into u and A'^i is approximately one. Therefore this 
gives the leading contribution to the process we are interested in. In principle, our 
calculation can be seen as a part of the longer decay chain, where the decays 
into two leptons and the LSP. This larger chain could be needed anyway, since 
outgoing muons are easy to detect which could reduce background. Here we do 
not consider the background from SM or other SUSY processes and the smearing 
of the signal from limited detector resolution. 

Of course there is no mixing between the processes with different incoming 
quark flavors. All contributions are simply added. But for each flavor interfer- 
ence terms in principle occur, after exchanging the b- and 6-jets. They are simply 
evaluated by Madgraph, by calculating the matrix element numerically and then 
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squaring it. We tried to compute the diagrams, including interference terms in this 
way. Unfortunately this integration is quite unstable and converges very slowly. 
This can be understood by having a look at fig. 16.11 When the momentum of 
the 6-jet from the decay of the quasi-on-shell gluino to the sbottom is calculated, 
it is found to be much softer than the 6-jet from the decay of the sbottom to 
the NLSP, as later shown in fig. 16.141 If these momenta are now exchanged for 
the calculation of the interference term, the gluino and sbottom can not both be 
on-shell. This makes our integration very inefficient because the Breit-Wigner 
mapping effectively only evaluates points at the pole. Since the interference terms 
are not increased by the Breit-Wigner propagator, they are suppressed in com- 
parison to the squared matrix elements without the exchange of the b- and 6-jet. 
This Breit-Wigner suppressed contribution to the angular distribution is expected 
to have the same shape as the pure squared matrix elements without interference 
terms and can therefore be left out as a NLO effect. As a consequence we calculate 
the diagrams with exchanged b- and 6-jets separately. We then find that then the 
integration is stable and converges much faster. 

Another problem occurs from the contributions deriving from different inter- 
mediate squarks. Of course, there is interference between the diagrams with the 
same final and initial states and including qi or in the intermediate states. But 
it is not trivial to calculate interference of such diagrams, since the phase space 
has poles from the Breit-Wigner propagators at different energies. This problem 
can not be solved in our calculation since the pole either is set to or rrig^. 
The phase space can only be generated for one of these poles at a time, e.g. . 
Then the resulting four momenta do not generate on-shell contributions to the 
other diagram with the pole at . The region of the phase space where many 
points are evaluated for a pole at niq^ is Breit-Wigner suppressed for the diagrams 
with the pole at rriq^. Therefore interference terms between these diagrams are 
again of NLO and therefore left out. Of course it is crucial for this approximation 
that the difference between rriq^ and m^^ is much larger than the width of the 
particles, as one can see in table ETJ The same argument is used for the sbottoms 
in the intermediate state and equally in the case of the UED scenario. 

If we want to exchange the A'^i and in the final state we would also have 
to neglect their interference term with the same argument since their masses are 
quite different. 

All these different contributions to the final state are added and filled into 
histograms. Altogether sixteen integrations of different phase spaces have to be 
evaluated, each summing over particles, antiparticles and the four light flavors 
and including the three different topologies as given in fig. 16.11 This results in a 
run time of several hours. For the case of UED we included the particle spectrum 
with quantum numbers and names of the particles into the file particles . dat of 
the Madgraph code. All Feynman rules needed for our process are included into 
the file iiiteractions.dat to generate the correct couplings. It is necessary to 
use the correct couplings, even if one normalizes all cross sections to one, since 
the mixing angle between singlets and doublets will be varied later. Since we 
want to exclude some interference terms as well as some topologies explicitly, as 
explained in the section 16.21 we can comment out the unwanted couplings in the 
interactions.dat file. 
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Figure 6.2: Diagrams contributing in LO to the 2^5 process in UED. Ah hght 
incoming quark flavors can not be distinguished and have to be summed over. 
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Figure 6.3: Diagrams are suppressed since the particles in t- and u-channels can 
not be on-shell. 

Using the quantum numbers of all particles as given in the Madgraph files, 
HELAS calculates the matrix element as the sum of combinations of helicity 
eigenstates of the external particles as we explained in section 15.21 Each SUSY 
contribution is calculated as the sum of 128 helicity combinations. In the case of 
the UED scenario there are 288 different helicity combinations. The number of 
helicity combinations is different, since the UED final state contains the massive 
gauge bosons ^3^1 and Bi, both having three helicity eigenstates. In SUSY, the 
fermions A^i and only have two helicity eigenstates. The boson Bi in UED is 
the KK-partners of the photon, corresponding to the SUSY LSP A^i, while the 
^3^1 corresponds to the NLSP Due to the larger number of helicity combina- 
tions the UED program has a much longer run-time because the matrix element 
routines are the slowest routines of the code. 



6.2 Neglected Topologies in the Final Decay Chain 

Using SUSY-Madgraph for the construction of the topologies and taking into ac- 
count all existing MSSM couplings, SUSY-Madgraph finds 460 different diagrams 
for the process denoted by the external particles 

u g ^ Ni ubb N2 ■ 

The same number of diagrams comes from 

dg^NidbbN2, 

where the quantum numbers of the quarks from the electroweak sector are differ- 
ent. Also antiquarks contribute to this process as we mentioned in the last section. 
In this section we shortly present the diagrams we neglect in our calculation. 

By using special examples from these Madgraph diagrams, we want to explain 
that diagrams different from those in fig. 16.11 can be neglected. Due to the Breit- 
Wigner propagators, s-channel-like decays give the dominant contribution to the 
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Figure 6.4: This diagram is suppressed since ul and g can not be on-sheh at the 
same time. 

resulting cross section. The diagrams presented in fig. 16.11 involve three Breit- 
Wigner propagators, all being on-shell in our calculation. 

As we already stated for the topologies of 2 ^ 3 and 2^4 processes, particles 
in t- and u-channels can not be on the mass shell. Therefore the topologies from 
fig- 16.31 are suppressed. There are a lot more combinations of particles which 
can be included into these topologies. Additional t-channel topologies also exist, 
including one or two 1 — > 2 particle decays on the right side of the t-channel. 

Another kind of topology that was neglected is given in fig. 16.41 Considering 
the mass spectrum of the parameter point SPS la, one finds that the gluino is 
heavier than all squarks. Therefore either the ml or the g can not be on-shell in 
fig. 16.41 Although only s-channel decays are involved, this diagram is suppressed 
by a missing on-shell Breit-Wigner. This constrains the validity of our calculations 
to mass spectra having the same mass hierachy. 

For a UED-like degenerate mass spectrum this topology could not be neglected. 
For the case of a nearly degenerate mass spectrum one also faces the problem that 
outgoing 6-jets are very soft and therefore hard to see at LHC experiments. 

Due to the gluons and quarks in the initial and final state there are different 
diagrams including gluon propagators. One of them is given on the left side 
in fig. 16.51 Since the gluon has no finite width included in its propagator, the 
diagram is suppressed. This holds for a large number of diagrams with internal 
gluons and only two additional Breit-Wigner propagators. Another example for 
that is given on the right side in fig. 16.51 Here even a four boson vertex is involved. 
This diagram would not occur in UED in the same way, since the 6i would be a 
fermionic partner of the bottom quark then. 

Of course, for our UED calculation we consider exactly the same topologies as 
in the case of SUSY. In principle, diagrams with higher excitations of KK-towers 
also exist. These are suppressed by the higher mass of the particles. 
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Figure 6.5: These diagrams are of NLO since there are only two on-shell Breit- 
Wigner propagators included. 

6.3 Results of the SUSY-UED Comparison 

To become more familiar with kinematics in decay chains we present some plots 
of invariant masses and momenta for UED and SUSY in this chapter. We also 
present the results for the angular distributions and the boost of the gluino/KK- 
gluon. Though it is in principle fixed by the masses, we use the mixing angle a^^\ 
introduced already in chapter O to test if the SUSY angular distributions can be 
reproduced by a UED decay chain with modified couplings. We therefore calculate 
the chain for 0° and 45° mixing angle of the KK-quark towers, beginning with the 
0° scenario. Both KK-towers in general can be expected not to mix due to the 
low mass of the SM 6-quark. All distributions are normalized to one in order to 
make both theories comparable. Since the error of our integration of the total 
cross section is below 1% there are no error bars presented in the histograms. 

6.3.1 A SUSY-UED Comparison for a^^) = 0° 

As we already stated in section [531 the invariant masses and transverse momenta 
and rapidities of intermediate particles served as a good check when building up 
our decay chain program stepwise. Since contributions nearly completely derive 
from on-shell momenta, these curves do not change for all steps of our decay chain 
programs. For example, the transverse momentum of the gluino stays unchanged, 
no matter what particles it decays to as long as the decay-vertex itself is not 
dependent on momenta. Since we did not present these distributions yet, we now 
want to present them for the 2 — > 5 process discussed in the last section. 

In fig. 16.61 we present the invariant mass distribution of the gluino and the 
KK-gluon for the SUSY and the UED decay chain. One finds two strongly peaked 
Breit-Wigner distribution exactly at the same value for the invariant mass. The 
invariant mass is given by 



^ij = \J{Pi + PjY ! (6-1) 

where pi and pj can be four vectors of two different particles. For the invariant 
mass of one single particle with momentum pi we use = yipip- 

For the invariant mass distribution of the sbottom h and KK-6-quark one finds 
similar distributions, shown in fig. 16.71 Depending on the mixing angle of the 
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Figure 6.6: Invariant mass distribution of the gluino and KK-gluon. 
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Figure 6.7: Invariant mass distribution of the 6-squark and KK-6-quark. 



two KK-towers in UED, either the first or second tower does not contribute to the 
2 — > 5 process when the mixing matrix is diagonal. Here we assume a mixing angle 
of a*^^^ = 0° between the two KK-towers and obtain a peak for the KK-quark 61^2 at 
516.91 GeV. This is also reflected by the Feynman rules in appendix lA. 4.31 While 
the coupling strength of one tower is at the maximum, the coupling to ^3^1 for the 
second tower is equal to zero. In the case of SUSY both sbottoms contribute to 
the process and therefore a second peak for sbottom 62 at 546.24 GeV is present 
in addition to the sbottom bi peak at 516.91 GeV. 

The invariant mass of the u,d,c and s-squark/KK-quark is given in fig. 16.8^ 
showing that ul does not decay to A''i with a large branching ratio. The two 
masses of qi^i and qi^2 in UED are chosen such, that the larger peak for UED 
is to be found at the same place as it is the case for SUSY. Also in the case of 
q;(o) = 0° mixing angle, there are two peaks for the UED scenario, since none of 
the couplings becomes zero. 

In fig. 16.91 we show the differential cross section for the transverse momentum 
of the gluino, of the light jet in fig. I6.1UI and of the non-measurable neutralino 
in fig. 16.111 Of course, neither the transverse momentum nor the invariant mass 
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Figure 6.8: Invariant mass distribution of the u,d,c and s-squark and KK-quark. 
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Figure 6.9: Transverse momentum distribution of the gluino and the KK-gluon. 
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Figure 6.11: Transverse momentum distribution section of A^i and Bi. 




Figure 6.12: Rapidity distribution of the gluino and the KK-gluon. 
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Figure 6.13: Distribution of the boost parameter /3 of the gluino and KK-gluon. 
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distributions of the intermediate gluinos/KK-gluons and squarks/KK-quarks can 
be measured directly in experiment. While the transverse momentum of gluino and 
KK-gluon are quite equal, there are differences between the transverse momentum 
distributions of the light quarks in both scenarios and between the transverse 
momentum distributions of A^i and Bi. These differences are due to the different 
rapidity distributions of the u,d,c and s-squark in SUSY and UED. This can 
already be observed for our 2^2 process in chapter [H The shape of the rapidity 
distribution close to y = differs for SUSY and UED and stays unchanged also for 
the 2 — > 5 process. The Ni, which is the LSP and the Bi, the lightest Kaluza-Klein 
particle (LKP), can not be observed directly. 

From the rapidity distribution in fig. 16.121 we also do not see significant differ- 
ences between gluino and KK-gluon. As we expect, these distributions are nearly 
equal to those of the 2 — > 2 process. 

As a result the distributions of the boost parameter /3 for the gluino and the 
KK-gluon, given in fig. 16.131 are quite similar. Here /? is given by 



One finds that the boost distribution for large values of /5 is slightly higher for the 
gluinos than for the KK gluons. Our plot of the (3 distribution qualitatively agrees 
with the one given in [49(]. There gluon-gluon collisions are taken into account 
additionally, while we only considered the process of quark-gluon collision. In 
the difference of the boosts seems to be more significant which could be due to the 
acceptance cuts assumed there. In our analysis we do not take into account any 
cuts or smearing of the signal due to finite resolution of the detector. In principle, 
cuts on the 6-jet momenta are necessary since the detector can not distinguish 
between both jets if they are very close to each other. One usually uses transverse 
momentum cuts which can have an impact on the shape of the curves. 

The 6-jet coupling directly to the gluino/KK-gluon, is usually called the near 
6-jet while the second 6-jet, emerging from the decay of the 6-squark/KK-6-quark 
is called the far 6-jet. They have, due to the mass hierachy in the decay chain, very 
different transverse momenta. Because of the invariant mass differences between 
gluino and sbottom and between sbottom and the NLSP, the second 6-jet is much 
harder than the first one. This results in very different transverse momentum 
distributions for the near and the far jet, as shown in fig. l6.14l for the UED scenario. 
A comparison of the near and far jet transverse momentum for SUSY and UED is 
shown in fig. 16.151 and fig. 16.161 One finds that the peaks are only slightly shifted. 

Due to different spins of the gluino and KK-gluon and the 6-squarks and KK- 
6-quark, we expect to find significant differences in the angle 6fji between these 
two 6-jets. Since in our program we boost all momenta to the lab frame, we show 
the angle between the two 6-jets in the lab frame in fig. 16.171 This angle Of^i is 
not Lorentz invariant under boosts along the beam axis. But this is true for the 
azimuth angle 0^^* between the two jets, i.e. the angle in the plane orthogonal to 
the beam axis. The azimuth angle is plotted in fig. 16.181 For the calculation of 
6fjl we used 




(6.2) 




with 0° < < 180' 



lO 



(6.3) 
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Figure 6.14: Transverse momentum distribution of the near and far 6-jets in UED. 
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Figure 6.15: Transverse momentum distribution of the outgoing near 6-jet. 
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Figure 6.17: Distribution of the angle between the two bottom jets in the lab 
frame. 



The same relation is used for 9^^^ with the third component of both vectors equal 
to zero, which is the direction of our beam axis. 

In both plots we find significant differences between the two models. It seems 
that a large part of this effect is not due to the different boosts which are quite 
close for both models. The origin of this difference in the angular distribution is 
discussed in more detail in the following sections. The angular distribution given 
in 49(] can be distinguished equally well, again containing quark- gluo n and gluon- 
gluon collisions. The plot for the azimuth angle looks different in [49| which could 
again be due to the cuts on the ^-jet momenta. To quantify the difference between 
the angular distributions in SUSY and UED the following asymmetry seems to be 
a reasonable measure: 



A 



ct(A(/.,-, < 90°) - aiA^,i > 90°) 
a(A0,5<9O°) + c7(A0,5>9O°) 



(6.4) 



For SUSY we obtain A = 0.1443 while for UED we find A = 0.0903 in the case of 
mixing angle a^^^ = 0°. 

As in [49! we also present the distribution of the averaged rapidity of the 
outgoing b- and 6-jets, given by 



V 



Vb + Vi 



(6.5) 



While they are different in [49|], we do not find a significant difference between 
them in fig. I6.19[ 

As a last plot we present the differential cross section for the invariant mass of 
the two bottom jets in fig. 16.201 We already mentioned that we assume all masses 
to be already measured in this thesis. At the LHC, masses are extracted from 
invariant masses of outgoing particles. As discussed in references [1^ and (4o| . 
the edges of invariant mass distributions will be used to measure the masses of 
the particles. In the case of the gluino decaying into two 6-jets and the NSLP, the 
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Figure 6.18: Distribution of the azimuth angle between the two bottom jets. 
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Figure 6.19: Distribution of the average rapidity of the outgoing 6/6-jets. 
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Figure 6.20: Invariant mass distribution of the 66-jet pair. 
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edge is given by 




(6.6) 



and equally in the case of UED. Using the SPS la mass spectrum, we find a value 
of nif^i = 296.62 GeV for the right edge with an intermediate bi. For the left edge 
with an intermediate 62 one finds rrij^i = 247.93 GeV. This agrees quite well with 
the edges in fig. 16.201 

The small amount of events lying on the right side of the edges is due to off 
shell effects, i.e. when the gluino/KK-quark or the sbottom/KK-6-quark is off- 
shell the 6-jet momentum can be increased. We checked that the number of points 
on the right side of the edge lowers if the integration parameter n, the number 
of widths integrated over as given in eq. (j5.38p . is reduced. Of course, there is 
only one KK-6-quark contributing to the invariant mass distribution in the case 



Until now we always assumed that there is no mixing between the two KK- 
quark towers. In the following we present the results for the maximal mixing 
of a(i) = 45°. 

One finds significantly different angular distributions. Since the couplings of 
all squarks change, the distributions of pt and y of the 6-jets also change. Their 
transverse momentum distributions are given in fig. 16.211 and fig. 16.221 Especially 
the transverse momentum distribution of the near 6-jet exhibit a different shape. 
In fig. 16.231 and fig. 16. 241 we find the differential cross sections for the angle between 
both 6-jets in the lab frame and the azimuth angle as defined before. Both curves 
can not be discriminated as in the case of 0° mixing angle. Obviously the coupling 
structure seems to have a strong impact on the angular distributions. 

The distribution of the boost parameter f3 is found to be equal to the case 
of a^^^ = 0°. Other distributions like invariant mass of the gluino do not change 
either. The invariant mass of the squarks is somewhat special since there are two 
masses that can be chosen in two different ways. For our calculations we choose 
the masses of 61^1 and 61^2 such that the invariant mass distribution in the case of 
Q,(i) = 0° has its highest peak at the same place as it is the case in SUSY. This is 
obvious from fig. 16.71 The same choice applies to the invariant mass of the u,d,c 
and s squark/KK-quark. Its invariant mass for a^^^ = 45° is shown in fig. 16.251 

The invariant mass distribution of the 6-jets for q(^) = 45° also looks signif- 
icantly different. In fig. 16.261 one finds that the main contribution to the UED 
cross section derives from an intermediate 61^2- But due to mixing there is also a 
contribution from an intermediate bi^i as opposed to the case of a(i) = 0°. 



of a(0) =0°. 



6.3.2 A SUSY-UED Comparison for a^^) = 45° 
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Figure 6.21: Transverse momentum distribution of the outgoing near 6-jet. 
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Figure 6.22: Transverse momentum distribution of the outgoing far 6-jet. 
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Figure 6.23: Distribution of the angle between the two bottom jets in the 
frame. 
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gure 6.24: Distribution of the azimuth angle between the two bottom jets. 
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jure 6.25: Invariant mass distribution of the u,d,c and s-squark/KK-quark. 
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6.4 Influence from the Boost of the Gluino/KK-Gluon 

As mentioned in 0] the difference in the angular distributions in fig. 16.181 could be 
due to the boost and the helicity structure of the squark couplings. To find out 
how strong the influence from the boost really is, we map the UED boost onto the 
curve of the SUSY boost by multiplying the internal Vegas phase space weights 
with the ratio of the height of the histogram bins from SUSY and UED. The same 
procedure is applied to all distributions, multiplying them for each point in phase 
space with the appropriate factor belonging to the bin of the boost at that phase 
space point. This is done for each phase space point evaluated by Vegas. Thereby 
we obtain all other distributions, especially the angular distributions, assuming at 
the same time that the boosts for the gluino and the KK-gluon are equal. The 
boost distributions are then by construction given by the SUSY curve in fig. 16.131 

The distributions for the azimuth angle between both jets are given in fig. 16.41 
We find that the influence from the boost is negligible because the boost of the 
gluino and the KK-gluon are quite similar. Therefore we expect the difference in 
the angular distributions to be mainly determined by the coupling structure of 
the KK-quarks and the different helicity eigenstates of the outgoing A'^i and 
respectively Bi and ^3,1. Other kinematic observables do not change either, when 
the UED boost for a^^^ = 0^°^ is mapped onto the SUSY boost. 

The strong influence from the coupling structure can also be seen in section [6^ 
where a mixing angle of q*-^^ = 45° is used. The angular distributions are much 
closer in this case, though then other kinematic distributions, like the transverse 
momentum of the near 6-jets, significantly differ. 
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Figure 6.27: Distribution of the azimuth angle between the two bottom jets. 
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Conclusions 



In order to distinguish between different scenarios of new physics at the LHC, it is 
important to measure as many properties as possible of all new particles. In this 
thesis we concentrated on the determination of the spin of an intermediate particle 
by using a decay chain. The gluino is especially interesting for SUSY searches at 
the LHC, since it is a Majorana fermion which can produce like-sign dileptons 
in the final state. In UED the like-sign dilepton signature could occur as well 
from the decay of a bosonic KK-gluon. Therefore it is crucial to measure the spin 
of the decaying particle to differentiate between different beyond-Standard-Model 
scenarios. 

In the case of SUSY, we calculated the partonic cross sections of all 2 — > 2 
SUSY-QCD processes. The quark- gluon collision was analyzed numerically on the 
hadronic level. After calculating the partonic cross section for the production of a 
KK-quark-gluon pair in UED, a comparison to the corresponding SUSY process 
was performed at the production threshold. It was found that the UED 2^2 cross 
section for quark-gluon collision increases much faster at the threshold, compared 
to the one in SUSY. In a comparison of the squark and KK-quark rapidities a 
slight difference between the curves could be identified. 

After explaining general kinematics in decay chains, we demonstrated how to 
test a decay chain program by using invariant masses, transverse momenta and 
rapidity distributions of intermediate particles. The shape of the curves of an 
intermediate particle in a chain does not change if the chain is extended by an ad- 
ditional decay. This is due to our integration procedure, only considering momenta 
close to the Breit-Wigner pole. Within this argumentation the properties of the 
Breit-Wigner propagator are crucial. Since a Breit-Wigner propagator effectively 
lowers the order in the coupling constant, some topologies can be neglected as they 
are effectively of higher order. Our program was compared to the narrow width 
approximation, i.e. only considering the on-shell contributions by multiplication 
of the particle production by the branching fractions. 

Hadronic cross sections were compared for the case of a well-known decay chain 
involving the SUSY gluino or a UED KK-gluon and assuming equal masses for 
the corresponding particles in both scenarios. In a comparison of this decay chain 
for SUSY and UED we could show that significant differences in the distributions 
of the azimuth angle between the outgoing 6-jets occur. As a consequence, we 
showed that it is possible to distinguish a SUSY gluino and a UED KK-gluon by 
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using angular correlations and assuming a naturally mass suppressed UED mixing 
angle of a^^) = 0°. 

By mapping the boosts of UED and SUSY onto each other, a comparison 
of the distributions of the azimuth angle between the 6-jets and other kinematic 
distributions is performed while effects from the difference in the gluino and KK- 
gluon boost are effectively eliminated. We found that the difference in the boost 
distributions is too small and its effect is not large enough to cause significant 
differences in the azimuth angle between the outgoing 6-jcts. 

For the outgoing light jet and the lightest supersymmetric particle, respectively 
the lightest Kaluza-Klein particle, we found differences in the kinematic distribu- 
tions. These derive from differences in the rapidity distribution of the squark and 
KK-quark in the final state of the 2^2 process. As we could show, this is most 
probably due to different angular distributions that can already be found on the 
partonic level in the 2 — > 2 process. 

It seems to be possible to make the azimuth angle distributions from SUSY and 
UED more similar by changing the mixing angle between the KK-quark towers 
for singlets and doublets. In the case of a^^^ = 45° the transverse momentum 
distribution of the near 6-jet shows significant differences for SUSY and UED, while 
the azimuth angles between the two outgoing 6-jets are nearly equal. Therefore it 
is most likely that the coupling structure and additional helicity states in the final 
state of the UED decay chain cause the differences in the angular distributions. 
The boost of the gluino and the KK-gluon remains nearly without any measurable 
effect. 

Though we could show that distinguishing the chains in a SUSY and UED 
scenario by kinematic observables is indeed possible, a more complete treatment, 
including off-shell effects, background effects and simulation of the finite resolution 
of the detector, would be needed to obtain a measurable result. Including next- 
to-leading order effects would improve the accuracy of the calculation. Especially 
QCD corrections can have a considerable size. 
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A.l Color Factors in Q CD-couplings 

Quantum Chromodynamics is a non-abelian gauge theory, i.e. there are symmetry 
transformations 

^ ^ e*"'*'V (A.l) 

which leave the theory unchanged. The are called symmetry generators. The 
symmetry in the color space of QCD has an underlying SU(A^) symmetry with 
= 3. Two of their representations, occurring in the QCD Feynman rules, are 
the fundamental and the adjoint representation. Quark fields transform under 
the former while gluon fields transform under the latter representation. When 
calculating QCD diagrams, sums over colors and flavors, i.e. traces over the sym- 
metry generators, appear in the formulae. We give the formulae which were used 
in the calculations of the partonic cross sections in chapter [H The generators of 
the fundamental representation are denoted as t". The appropriate Lie-algebra is 
given by 

[e, t^] = et^ - th" = ir^H^ , (A.2) 

where f"'^'^ are the structure constants of the Lie-algebra. In the following formu- 
lae, the lower indices of are always implicit. 

trifthH^) = tr (^{C2{r) - ^^2(0)^ tH"^ = 

/ Mbc ■ 

tr{tHHh^) = 3 {C2{r) f = y (A.3) 

jacdjbcd ^ C2iG)6"-^ 

tX = C2(r) • 1 

with 1 being the d{r) x d{r) unit matrix and C2{r) and C2{G) being the quadratic 
Casimir operators of both representations. The structure constant f"'^'^ is totally 
antisymmetric while d"'^'^ is totally symmetric. 
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The Casimir operators for the fundamental representation are given by 



CiN) = - 



C2{N) 



(A.4) 



2N 



The Casimir operators for the adjoint representation are given by 



C2{G) = C{G) = N . 



(A.5) 



The so cahed Jacobi identity gives 




(A.6) 



A. 2 An Explicit Calculation with External Gluons 

When more than one external gluon is involved in a calculation, one has to prop- 
erly take into account their non-physical longitudinal degrees of freedom. These 
longitudinal states come from the fact that Quantum Chromodynamics (QCD) 
is a non abelian gauge theory. They appear in the path integral quantization of 
QCD as non-physical fields usually called ghosts, spin zero fields with Fermi-Dirac 
statistics. In QED path integral quantization ghost fields also appear but decouple 
from the rest of the theory. 

Since calculating squared matrix elements by hand is a tedious work if one has 
to use the transverse polarization sum given in eq. (j4.4|) . it would be nice if one 
could simply use a shorter polarization sum equal to the one in QED, —g^u- In 
general the arbitrary vector has to drop out in the end since the result can 
not depend on this arbitrary, non-physical object. The result has to be gauge 
independent. Since the additional terms in the polarization sum given in eq. (j4.4p 
take out the longitudinal degrees of freedom, one has to do this by hand, if —gfj,u 
is used instead. 

In leading order matrix elements this is possible by using a trick called ghost 
subtraction. Here the non-physical polarizations of the gluon are explicitly sub- 
tracted from the matrix element by crossing out all terms proportional to momenta 
which have the Dirac index of the polarization vector and the same momentum. 
These are exactly those terms in the matrix element that drop out when eq. (j4.4p 
is used. It corresponds to using eq. ()4.6p already at the amplitude level before 
squaring the matrix element, i.e. a projection on physical degrees of freedom. The 
only exception is the case that all external particles are gluons. Then one still has 
to use eq. (|4.4p for one of them. 

Longitudinal degrees of freedom are taken out in order to make the (on-shell)- 
Ward-identities from QED also hold for QCD. liM{k) = e^{k)M^'ik) is the ampli- 
tude for a given QED process with an external photon e^, the amplitude vanishes 
if is replaced by /c^, i.e. it stays unchanged under the gauge transformation 
^fi = + k^. The so-called Ward identity is then denoted by 



In the case of QCD the equivalent relations are called Slavnov- Taylor identities. 

Following an argument in ^j, we want to explain how ghost subtraction ex- 
plicitly works. Therefore we check the Ward identity explicitly for the example 



k^Mf'ik) = . 



(A.7) 
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of gluino-gluino production from gluon-gluon collision. In this case the physical 
incoming gluons must have transverse polarization since they are on-shell. The 
matrix elements on Born level for the three contributing diagrams in fig. 14.11 e) , 
are given by 

^-^l.diagrame/.(^l)e^(^2) = 

= {-gs)r"'[g^''{kr - k^Y' + g"""' {k2 + ki + k2r + gP'f^{-ki -k2- kiY] 

(fci + k2y 



^-^2.diagran/M(^l)er.(^2) = 

— mi + i€ 



iM'i''^:^„ e^{ki)euik2) 



3. diagram M 



pr2 _ j^2 _j_ 



l''{-gs)r'v(j,i)e,{k2)e^{k,) 



with 



p = ki — pi and P = ki — p2 ■ 
The last two diagrams sum to 

iM^''Mki)eAk2) = u{p2) Wi-gs)!''" '^^ t '^'^2'' ri-gs)r'+ 

+ 7^/^'^(-gs) ""i^"' (-gs)r-^'^j v{pi)e,{k^)e,{k2) . 

Using eq. (j4.6p and for the gluinos 

(p — in)u{p) = and v{p) {—p — m) = , (A. 8) 

and replacing ey{k2) by k2u one obtains 

iM>i'^^e^{ki)k2, = -igsu{P2)[f'"T' - f'n'jre^.iki) . 

This term needs to be canceled by the contribution of the first term since the 
on-shell Ward identity has to hold. For the first term one finds 

1 



*^^rdiag.^M(fci)fc2. = -^g'r'''f'' \^^^^^,2 [a'^kl - g^^kl - k{e, + kP,k>^] 



(ki + k2)^ 
X u{p2)Yv{pi)f-i^{ki) 



(A.9) 



with ^3 = (fci + k2)- For the third term in brackets we now assume transversality 
of the external gluons. Therefore it is equal to zero when it is contacted with the 
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polarization vector. The first term in brackets disappears for on-sliell incoming 
gluons and the last term vanishes when it is contracted with fermionic currents, 
due to eq. ()A.8p . Therefore one obtains the result 

iMl':iiagrMh)k2u = ^5.' (" l)n(p2)7^^(Pl (^1 , 

which with the Jacobi identity from eq. ()A.6|) exactly cancels the contribution 
from the last two diagrams. 

Therefore the on-shell Ward identity is fulfilled and shows that it is indeed 
possible to treat the incoming particles as exclusively transverse polarized gluons 
by using ghost subtraction. Now —gi_iu may be used instead of the full polarization 
sum, since longitudinal degrees of freedom are no longer included in the matrix 
element. Using the Feynman rules for the ghosts given in one finds that the 
third term in brackets in eq. (lA.Op exactly cancels against the ghost fields if they 
are taken into account explicitly and cancellation is not performed by hand. 

When we used this method in section 14.21 we also checked gauge invariance 
of the result by using the full polarization sum from eq. ()4.4p and calculating the 
result in a general gauge without making use of the transversality condition, given 
in eq. (fOj) . 
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A. 3 SUSY-QCD Fey nman- Rules 



— ml + ie 



— + ie 




a, a p — > P,b 




a, /i 



-i9s{t'')ji{pi-PjT, 






±^(t'^),,(l±75W, 



±^(t'^),,(l±75W, 



For the last four graphs, the upper sign gives the Feynman rules for a left handed 
squark and the lower sign for a right handed squark. 



88 



Appendix 



A. 4 UED Fey nman- Rules 

A. 4.1 UED-Interactions with Gluons and Quarks 



a,i q b, j 



i{f + rnq)ba5ij a,i q„.^ b,J ^ + Mn)baSr 

— ml + ie ' 



P 





Qn,l 



cb ■ 



/ qn.:2 

\ 9? 1, 2 



cb : 



)WVW ^■2n*=ig_^>^^^(^t'^)^^^ 



' Qn,2 
Qn,2 



QnA 



;::/vvvw ^2n* = ")/vvvw ^"*= -i5-=7'^(t")^i, 



'V2 



V2 



/" Qn,2 



/ Q2n,l 



,.•■•■"9)71,2 



9771,1 

9ri.l 



V2 



■ 9777,1 



-ig—r{e)^,, 



Indices for all vertices are equal to those at the quark-quark-gluon vertex. The 
coupling is given by o = For better readability we use different lines for the 

V 7r/i 

heavy KK-partners of the gluons 5* and of the quarks The index z = 1,2 
denotes the heavy quark tower, while n denotes the excitation level within the 
tower. 
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A. 4. 2 Purely Gluonic UED-Interactions 



H, a gjj* i^, b 

/i, a Qj^^ ly, b 
WWWW 

k, fi, a 



+ (1 - A) 



p^ + ie 

Ml 



p^-Ml + i MnTl 



p^ + ie 




9 [{p - ky g^'^ + {q- pT g""' + {k- qT g'^] , 



A;, /i, a 



p, v,b 




V2 



+ f^"''^ f^'"^ {9^^ 9^" — 9^'^ 9^'^) 
+ f^""^ f^'^'^ {9^^ 9^" — 9^" 9^^)] 



9 
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gn* 



gn* 




93n* 



-I — q 
2 ^ 




g 




-I - q 
2 ^ 




-I - q 
2 ^ 



Z±m±n| 



All momenta are defined incoming. The indices belonging to the four particle 
interactions are the same as for the four gluon vertex. The coupling is given by 
. For better readability we use curved lines for the heavy KK-partners of 



g 
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the gluons (7*. The three and four particle vertices have to be multiplied by the 
tensor structure term in brackets, belonging to the three and four gluon interac- 
tions. 
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A. 4. 3 Selected Feynman Rules from the Electroweak UED Sector 




The Lagrangian of the electroweak sector is given in chapter 3. The angle a^^^ 
denotes the mixing angle of singlets and doublets of the first excitation level. The 
particles qi^n here denote the mass eigenstates, derived from the eigenstates of 
weak interaction by mixing of the two towers. Additional electroweak Feynman 



rules can be found, e.g. , in [26(] 
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